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Abstract 

In this paper, a novel concept called a uniquely factorable constellation pair (UFCP) is proposed 
for the systematic design of a noncoherent full diversity collaborative unitary space-time block code by 
normalizing two Alamouti codes for a wireless communication system having two transmitter antennas 
and a single receiver antenna. It is proved that such a unitary UFCP code assures the unique identification 
of both channel coefficients and transmitted signals in a noise-free case as well as full diversity for the 
noncoherent maximum likelihood (ML) receiver in a noise case. To further improve error performance, 
an optimal unitary UFCP code is designed by appropriately and uniquely factorizing a pair of energy- 
efficient cross quadrature amplitude modulation (QAM) constellations to maximize the coding gain 
subject to a transmission bit rate constraint. After a deep investigation of the fractional coding gain 
function, a technical approach developed in this paper to maximizing the coding gain is to carefully 
design an energy scale to compress the first three largest energy points in the corner of the QAM 
constellations in the denominator of the objective as well as carefully design a constellation triple 
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forming two UFCPs, with one collaborating with the other two so as to make the accumulated minimum 
Euclidean distance along the two transmitter antennas in the numerator of the objective as large as 
possible and at the same time, to avoid as many comer points of the QAM constellations with the 
largest energy as possible to achieve the minimum of the numerator. In other words, the optimal coding 
gain is attained by intelligent constellations collaboration and efficient energy compression. Computer 
simulations demonstrate that error performance of the optimal unitary UFCP code presented in this 
paper outperforms those of the differential code and the SNR-efficient training code, which is the best 
code in current literatures for the system. 



Index Terms 

Cross QAM constellations, constellations collaboration, coding gain, energy scale, full diversity, 
noncoherent ML receiver, uniquely factorizable constellation pair and unitary space-time block code. 



I. Introduction 

At present, the technology inteUigently combining multiple antennas |[T|-||3| with space- 
time block coding p]|-p5| has been well developed to improve the spectral efficiency of a 
coherent wireless communication system. Lately, simple space-time block code designs with low 



complexity decoding [24|-[39| have attracted much attention. In this paper, we are specifically 
interested in a flat fading wireless communication system with two transmitter antennas and a 
single receiver antenna. This system is often encountered in mobile down-link communications 
for which the mobile receiver may not be able to deploy multiple antennas. For such a system, 
if the exact knowledge of the channel coefficients is available at the receiver, the orthogonal 



Alamouti [26| space-time block code is particularly appealing, since it enables the coherent ML 
receiver to extract full diversity not only with linear processing complexity, but with information 



losslessness as well [40|. Unfortunately, perfect channel state information at the receiver, in 
practice, is not easily obtainable. If the channel changes slowly, then, the transmitter may 
have sufficiently long channel coherence time and send training signals enabling the channel 
coefficients to be estimated accurately. However, the fading coefficients in mobile wireless 
communications may vary rapidly and the coherence time may be too short to allow reliable 
estimation of the coefficients. Therefore, the time cost on sending training signals cannot be 



ignored since more training signals need sending for the accurate estimation of the channel |41 1- 



[43 1 . In this paper, we consider the communication scenario where channel fading changes 
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very promptly, assuming that the channel gains are completely unknown at both the transmitter 
and the receiver, but remain unchanged within four transmission time slots, after which they 
change to new independent values that are fixed for next four time slots, and so on. Such fast 
varying channel with flat fading for a single transmitter and single receiver antenna was first 



considered for determining the capacity-achieving input distribution [44|-[48|. In order to make 
communication as reliable as possible under this severe environment and avoid sending the 
training signals for estimation of the channel, using differential space-time block coding p9|- 



|59| is one of the possible solutions. Unfortunately, this approach results in an approximate loss 
of 3dB in performance compared to coherent detection. Recently, some techniques of blind signal 
processing such as the subspace method based on the second order statistics have been utilized 



to blindly identify the space-time block coded channel [60|-[64|. However, phase ambiguity 
incurs the channel not being able to be identified uniquely, even in the noise-free case. In 
addition, even if there were no phase ambiguity, the subspace method could not be successfully 
applicable to our case, since the 4 coherence time slots are too short to allow the second-order 
statistics to be estimated accurately. Therefore, in order to attain a more satisfactory solution. 



noncoherent space-time block coding techniques [ [49| , [|65||-[|69| have been developed. It has been 
proved that either at high signal-to-noise ratio (SNR) or for long coherence time, the unitary 



code is optimal [43 1, [49 1, [70 1, [71 1. Hence, most of the noncoherent space-time block code 



designs have been primarily concentrated on unitary designs [49|, [58[, [65|-[69|, [72[. The 



Cay ley [ [54| , [ |58| transform and the exponential transform [ [69| are now two well-established 
transforms which convert a linear dispersion code and a linear space-time block code into unitary 
codes. The exponential transform [ [69| requires that the number of the receiver antennas is greater 
than or equal to that of the transmitter antennas. In general, it cannot assure full diversity for 
the noncoherent ML receiver. The unitary design using the Cayley transform aimed mainly 
at differential modulation and a differential receiver. Recently, the original non-full diversity 



design [54|, [58 1 based on the Cayley transform has been improved into a full diversity design 
utilizing division algebra and algebraic number theory [ |73| . More recently, the systematic design 
of nocoherent unitary space-time block codes with full diversity and a high transmission rate for 
an arbitrary number of the transmitter antennas and the receiver antennas has been established 
by using a pair of coprime PSK constellations and the QR decomposition [)74|. Particularly for 
the system with two transmitter antennas and a single receiver antenna, the phase ambiguity and 
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full diversity issue for the noncoherent Alamouti space-time block code has been completely 
resolved [|75|, [|76l. 

However, the PSK constellation is not as energy-efficient as the QAM constellation. Therefore, 
our primary target in this paper is to design a full diversity unitary space-time block code for the 
system by using the two Alamouti codes and the energy-efficient cross QAM constellations such 
that the noncoherent coding gain is maximized subject a transmission bit rate constraint. Despite 
the fact that recent research on coherent MIMO communications has told us that the Alamouti 
code enables coherent full diversity for any constellations and with any receivers, this is no 
longer true for noncoherent communications, even if the commonly-used QAM constellations are 
transmitted and even if the noncoherent ML receiver is employed, since the likelihood function 



is invariant under certain rotation of some QAM constellation points [64|, [77|. Hence, signals 
must be carefully designed to combat against fading. In the noncoherent wireless communication 
scenario, the unknown of the fading channel at both the transmitter and the receiver requires that 
the transmitted signals emitting from different time slots must be more correlated than in the 
coherent environment so that reliable communications with noncoherent full diversity are made 
possible under a maximum allowable transmission date rate. 

All the aforementioned factors greatly motivate us to proposing a novel concept, a uniquely 
factorable constellation pair (UFCP), for the systematic design of an optimal unitary constellation 
for the system. The main idea of the UFCP design essentially comprises the following two major 
steps. 

1) Intelligent constellations collaboration: From a pair of the energy-efficient cross QAM 
constellations, a constellation triple constituting two UFCPs will be carefully designed, 
with the one collaboratively shared with the other two through the two transmitter antennas, 
so that the minimum of the numerator in the factional objection function is made as large 
as possible and at the same time, the largest energy points of the QAM constellations as 
many as possible are avoided to reach the minimum of the numerator. 

2) Efficient energy compression: An energy scale will be carefully designed to compress 
the first three largest energy points of the QAM constellations in the denominator of the 
objective. 

Notation: Most notations used throughout this paper are standard: column vectors and matrices 
are boldface lowercase and uppercase letters, respectively; the matrix transpose, the complex 
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conjugate, the Hermitian are denoted by (■)^, (■)*, (■)^, respectively; In denotes the N x N 
identity matrix; Notation Tr(M) denotes the trace of an K x K matrix M, i.e., Tr(M) = 
J2f=i ^ii^ whereas notation det(M) denotes the determinant of M; Notation ||M||f denotes the 
Frobenius norm of M; Notation f{x) = o{g{x)) denotes Umj._>.oo = 0; $ denotes an empty 
set. 

II. Channel Model and Unitary Space-time Block Coding 

In this section, we first briefly review the channel model in which we are interested in this 
paper. Then, we propose our transmission scheme and unitary code structure. 

A. Channel Model 

Let us consider a wireless communication system having two transmitter antennas and a 
single receiver antenna. The transmitted symbols from the two transmitter antennas arrive at the 
receiver via two different channels hi and h-z. Then, the discrete baseband received signal r can 
be represented as 

r = hiSi + h2S2 + i- (1) 

Throughout this paper we assume hi and /i2 are samples of independent circularly symmetric 
zero-mean complex white Gaussian random variables with unit variances and remain constant 
for the first 4 time slots, after which they change to new independent values that are fixed for the 
next 4 time slots, and so on, the explanation for which will be given in the ensuing subsection, si 
and S2 are two corresponding transmitted symbols from these two antennas, and ^ is a circularly 
symmetric complex Gaussian noise with zero mean and variance a^. 

B. Unitary Space-Time Block Codes 

Let A, Bi and B2 be three constellations to be designed. Then, our unitary space-time block 
code for the channel model ([T]) is basically generated by normalizing two Alamouti codes 
and is described as follows: First, randomly, independently and equally likely choose three 
symbols a E A,bi E Bi and 62 G B2 and then, transmit their normalized version from the two 
transmitter antennas within four time slots. During the first time slot, we transmit the signals 
Si = a/ i/|aP + + and S2 = in ([T]) at the same time from the respective Antennas 
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1 and 2. During the second time slot, the signals si = and S2 = a*/^y\a\'^ + + are 
simultaneously transmitted from Antennas 1 and 2, respectively. Collecting these two received 
signals yields 



ri 
r2 



1 



a 



VW+W+W \ a* 



or equivalently, 



hi 



6 

6 



(2a) 



(2b) 



V|aP + |6i|2 + |62P' 

where = (ri,r2)^,A = diag(a, a*),h = {hi,h2)^ and = (^1,^2)"^- In the rest two time 
slots, the second and third symbols are transmitted using the Alamouti coding scheme, i.e.. 



^^3 



|a|2 + |6i|2 + |62P 



hi 62 
-hi hi 



hi 
h2 



6 

^4 



or equivalently. 



by 



(2c) 



(2d) 



where n = (rs^r^y,^^ = (6,^4) and 



B 



hi &2 
-hi hi 



hi e Bi, 62 G B2. 



By stacking the above four received signals ( [2b| ) and ( [2d| ), the relationship between the transmitted 
and received signals within the four time slots can be represented in a more compact matrix 
form as 



r = Uh + ^, 
where r = iri,r2,rs,nf = (^1,6,6,6)'^ and 



(3a) 



U 



A 



fa \ 

a* 

hi &2 

V -b*2 bl ) 



(3b) 



for a E A, hi E i3i,&2 G B2. We would like to make the following comments on the unitary 
space-time block coded channel model (|3]): 
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1) Why do we need four time slots! For a noncoherent MIMO system with M transmitter 



antennas, it has been proved [71 1 that a necessary condition for a space-time block code to 
enable full diversity for the noncoherent ML receiver is coherent time T > 2M. Particularly 
for M = 2, we should need at least 4 time slots. Actually, in this paper, we consider the 
shortest coherent time slots enabling the unique identification and full diversity. See more 
details in Section |Vl 

2) Why is the symbol rate of | per channel use reasonable! The answer to this question is 
mainly motivated from the following observation: For a noncoherent MIMO communi- 
cation system with M transmitter antennas and receiver antennas, Zheng and Tse [4] 
have proved that in a high SNR regime and for the Rayleigh-faded channel, the average 
channel capacity is given by 

C = M* (^1-^^ log SNR + 0(1), (4) 

where M* = min{M, A^, [|^J } and T is coherent time. This benchmark result tells us that 
the original noncoherent MIMO system can be asymptotically regarded as M* (l — ^) 
parallel spatial channels and thus, the number M* (l — 4^) is the total number of degrees 
of freedom to communication. The result also suggests us that the symbol rate of a space- 
time block code for the noncoherent MIMO channel should be M* (l — ^!-). Especially 
for the noncoherent system with T = 4, M = 2 and A^ = 1, the symbol rate should be |. 



C. Problem Statement 

To formally state our design problem, we make the following assumptions throughout this 
paper: 

1) The channel coefficients hi and /i2 are samples of independent circularly-symmetric white 
Gaussian random variables with zero mean and unit variances, and remain constant for 
the first 4 time slots, after which they change to new independent values that are fixed for 
the next 4 time slots, and so on. 

2) The elements of ^ are circularly-symmetric zero-mean complex Gaussian samples with 
covariance matrix a'^\i, 

3) During 4 observable time slots, the space-time block coding matrix U is transmitted 
with a, hi and 62 being independently and equally likely chosen from the respective 
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constellations A, Bi and B2- 
4) Channel state information is not available at either the transmitter or the receiver. 
Under the above assumptions, our primary purpose in this paper is to solve the following problem. 

Problem 1: Design the constellation triple A^Bi and B2 for the unitary space-time block 
coded channel ([3]) such that 

1) in the noise-free case, for any given nonzero received signal vector r 7^ 0, the equation 



reduced from (3a) 



r = Uh (5) 

with respect to the transmitted symbol variables a, hi and 62, and the channel vector h 
has a unique solution, and 
2) in the noisy environment, full diversity and the optimal coding gain are enabled for the 
noncoherent ML receiver. ■ 

III. Uniquely Factorable Constellation Pair 

In order to systematically design the constellation triple. A, Bi and B2 in Problem [T] in this 
section we propose a novel concept called Uniquely Factorable Constellation Pair (UFCP). 

A. UFCP 

Definition 1: A pair of constellations X and y is said to be a UFCP, which is denoted by 
y X, if there exist x,x E X and y,y E y such that xy = xy, then x = x, y = y. ■ 



In Section IV, we will see that it is this kind of the unique factorization of constellations that 
enables the unique identification of the channel and transmitted signals as well as full diversity. 
The following example provides us with a trivial UFCP. 

Example 1: For any set y, if we take X = {!}, then, X and 3^ form a UFCP. ■ 
Example [T] tells us that the constellation pair based on the training transmission scheme naturally 
forms a UFCP. In this paper, we are interested in the design of non-trival UFCPs each element 
of which is a complex integer. To do that, we need to develop a necessary condition which a 
UFCP must satisfy. 

Proposition 1: Let X and 3^ form a UFCP If > 2, then, 0^ X. ■ 
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Proof: Since |3^| > 2, there exist two elements ?/i,y2 G y such that yi ^ ?/2- Now, suppose 
that E X. Then, we would have x t/i = x ?/2 = 0, which contradicts with the assumption 
that X and y constitute a UFCP. This completes the proof of Proposition [T] □ 
From Proposition [T] and Definition [Tj we can immediately obtain the following proposition: 
Proposition 2: For a pair of given constellations X and y with each having finite size and 
^ Af, if a new constellation Z is defined as 

Z=\^z:z='^,xeX,yeyy 

then, such a pair of X and 3^ constitutes a UFCP if and only if 

\Z\ = \X\ X |3^|. (6) 

■ 

Proof: Define a map r from X x y to Z as 

T : {x,y) e X xy = z e Z. (7) 

X 

It can be verified that r is a one-to-one correspondence if and only if ^ ^ |, i.e., xy ^ xy 
for {x,y) 7^ ix,y), which is equivalent to saying that a pair of the constellations X and y 
constitutes a UFCP. On the other hand, since both X x y and Z have finite size, r is one-to- 
one correspondence if and only if \Z\ = \X x y\ = \X\ x \y\. This completes the proof of 
Proposition [2} □ 

Proposition [2] tells us that a UFCP y X can be constructed by factorizing a constellation Z 
in such a way that each fraction is unique. For notation simplicity, this kind of construction is 
specifically denoted hy Z = ^. In general, when Z is sizable, it is not easy to find a non-trival 
unique factorization. However, for some special constellations, we can utilize this factorization 
to systematically construct a UFCP. 

Definition 2: A constellation is said to be rotation-invariant with respect to e^^ if every element 
in the constellation multipled by e^^ still belongs to the constellation. ■ 

Example 2: The square QAM constellation is rotation-invariant with respect to e^^ with the 
angle 9 being equal to 7r/2, vr and 3n/2. ■ 
In general, we always have the following property, which can be verified directly by the definition 
and thus, whose poof is omitted. 
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Proposition 3: Let Z be rotation-invariant with respect to e^'^^'^, e^'^ , e^^'^^'^ and 0,r,rj ^ Z, 
where r is real. Then, the following statements are true: 

1) If we let 

X = {1,J}, 

y = {Z : Z = Zre+ j^im G 2, Zre, Zi^ > 0} U {z : Z = Zre + jZim G Z, Z^e, ^im < 0}, 

then, such a pair of X and y constitutes a UFCP and Z = ^. 

2) If we let 

X = {1,-1, J, -j}, 

y = {Z : Z = Zre+ j^im G Z, Z^e, ^im > 0}, 

then, such a pair of X and y forms another UFCP and Z = ^. 

m 

The following two examples show us how to obtain the UFCPs by factorizing the 16-QAM 
constellation by Proposition [3j 

Example 3: Let Z be the 16-QAM constellation. By Proposition [3| a UFCP 3^ ~ A" is obtained 
by factorizing Z: 

y = {3 + 3j, 3 + j, l + 3j,l+j, -l-j, -3-3j,-3-j,-l-3j}. 

■ 

Example 4: Again, applying an idea similar to Proposition |3] to the 16-QAM constellation 
yields another UFCP: 

X = {1,-1, J, -j}, 

y = {3 + 3j,3 + j,l + 3j,l+j}. 

■ 

In addition, it is noticed that for a given UFCP y X, when the symbol x is fixed, all the 
fractions of the form z = -,y E y can be chosen only from a certain subset of Z, which will 



be employed in the design of the constellation triple A, Bi and B2 in Section IV Therefore, we 
particularly give a formal definition as follows: 
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Definition 3: Given a UFCP y X and a fixed x E X, a set generated from x, denoted by 

Z, = {z:z=y-,yey} (8) 

X 

is called a Group-x. ■ 
The groups have some interesting properties: 

Proposition 4: Let Z = ^. Then, the following three statements are true: 

1) Non-intersection: For any xi,X2 E X,xi^ X2, there is no intersection between Group-Xi 
and Group-X2, i.e., 

z,^ n z,^ = (9) 

2) Decomposition: The union of all the groups is equal to the original constellation Z, i.e., 

U,^x 2. = Z. (10) 

3) The number of groups is equal to and \Zx\ = \y\ for any x E X. ■ 
Proof: Statement 2) is derived directly from the definition of the Group. Here, we only examine 
Statements 1) and 3). For any xi, X2 E X and xi ^ X2, suppose that there exists some z belonging 
to Zx^ n Zx2- Then, z E Z^^ and z E Z^^ and thus, there exist y\,y2 E y such that = 
which implies that X\y2 = X2yi. Since X and y form a UFCP, we have xi = X2 and yi = y2, 
which contradicts with the assumption that xi ^ X2- Hence, Z.j,^ fl Z.^.^ = i.e., Statement 1) 
is true. 

Using the Statements 1) and 2), we obtain 



By the definition of the group, we have \Za:\ = \y\. Substituting this into ( [TT| ) yields |A'||3^| 



\Z\. This completes the proof of Statement 3) and thus, of Proposition |4j □ 

B. Unique factorizations of the modified cross QAM constellations 

In spite of the fact that Proposition [2] tells us that a UFCP y X can be constructed by 
factorizing a constellation Z, in general, the UFCP y X so derived from the given Z is 
not unique. In other words, the same constellation Z can generate two different UFCPs. For 
instance, consider an example below: 
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Example 5: Let Z be the 16-QAM constellation. If we let 

y = {3 + 3j,3-j,-l + 3j,-l-j}, 

then, it can be verified that such a pair of constellations X and 3^ also forms a UFCP, which is 
different from the one given in Example |4j ■ 
Now, a natural question is: for a given constellation Z, which pair of constellations X and 
y generated by Z is better? A general answer to this question is hard to be given. However, 
in this paper, we are interested in the UFCPs which are derived from the cross 2^-ary QAM 



constellation [78 1 and every element of which is still a complex integer. This implicitly requires 
that each element of X must be of the unit-norm, i.e., ±1 or ±j. Specifically, we focus on such 
a UFCP generated from the cross QAM constellation that the minimum distance of y is as large 
as possible. Since the conventional 8-QAM constellation does not satisfy the rotation-invariant 
property, we need to modify it into a new 8-QAM constellation so as to be rotation-invariant 
under e^^/"^, e^'^ and e-''^'"/^. For discussion self-containment, a formal definition of a 2^-ary 
modified cross QAM constellation is provided here. 

Definition 4: A modified 2^-ary QAM constellation Q is defined as follows: 

1) If K is even, Q is the standard square 2^-ary QAM constellation, i.e., 

Q = |(2m - 1) + {2n - l)j : -2^ + 1 < m,n < 2^|. 

2) If K = 3, Q is a new 8-QAM constellation modified from the conventional 8-ary QAM 
constellation, i.e., 

Q = {l + 3j, 1 + J,3 - J, 1 - J, -1 - 3j, -1 - J, -3 + J, -1 + j}. 

3) If K is an odd number exceeding 3, Q is the union of a horizontal rectangular QAM 
constellation and a vertical rectangular QAM constellation, i.e., 

Q = {{2m - 1) + {2n - l)j : -3 x 2^ + 1 < m < 3 x 2^, -2^ + l<n< 2^} 
[j{{2m - 1) + {2n - l)j : -2^ + 1 < m < 2^,-3 x 2^ + 1 < n < 3 x 2^}. 



September 21, 2011 



Submitted to IEEE Trans, on Information Theory 



XIA, ZHANG AND DUMITRESCU, ENERGY-EFFICIENT FULL DIVERSITY COLLABORATIVE UNITARY CODE DESIGNS 



13 



Proposition 5: Let Z be the given modified 2^-cross QAM constellation. Then, subject to 
X C {1, — 1, j, — j} with a fixed size greater than one, one solution to the following optimization 
problem: 



{'^opt, O^opt} = argmax min \yi - 1/2 1 (13) 



^=z 2/1 7^2/2 



is given as follows: 
1) If I A" I = 2, then, 

a) Fori^' = 3, 

b) Fori^ = 5, 



= {l + 3j,-l-3j,-l-j,l+i}. (14) 



= {-1 + 5j, 3 + 5j, -3 + 3j, 1 + 3j, 5 + 3j, -5 + j, -1 + j, 3 + j, 

-3 - i, 1 - i, 5 - i, -5 - 3j, -1 - 3j, 3 - 3i, -3 - 5j, 1 - 5j}. (15) 

c) For i^' > 4, the optimal y^^^^ is determined as follows: 

i) When K is even, 

3;^JJ = |(2f - 1 - 4m) + (2^ - 1 - 4n)j : < m,n < 2^ - l| 

|j|(2f -3-4m) + (2f -3-4n)j : < m, n < 2^ - l|. (16) 

ii) When K is an odd number exceeding 5, 





- 1 - 4m) + (2^ - 1 • 


Jf-3 K — b 

1 m=2 2 -l,n=3x2 2 -1 

- 4n)i [ 

J m=0,n=0 


U{(2^-i 


- 4m) + (3 X 2^ - 1 


m=3x2 2 -l,n=2 2 -1 

- 4n) j \ 

) m=0,n=0 


|j{(3x2^ 


-3- 4m) + (2^ -3 


m=2 2 -l,n=3x2 2 -1 

- 4n) j ^ 

J m=0,ra=0 


U{(2--3- 


-4m) + (3x 2^ -3- 


K-6 fsr-3 

>.m=3x2 2 -l,n=2 2 -1 

- 4n)j ^ 

J m,=0,n=0 


4 and > 3, then. 








^opt l-*-' -"-'J' 





(17) 
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a) For is: = 3, 



b) For if = 5, 



ySt = {i + 3j,-i-j}. 



yopt = {-l + 5j,3 + 5j,-5 + j, -l+j,3 + j,-5-3j,-l-3j,3-3j}. 

c) For K > 4, the optimal 3^Qp|^ can be determined as follows: 

i) When K is even, 

37^2) ^ |(4m - 2^ +3) + (2^ - 1 -4n)j : 0<m,n<2^-l|. (18) 

ii) When K is an odd number exceeding 5, 

( , K-3 , , K-i , ^ ' 

3^opt = (3 X 2^ - 1 - 4m) + (2^ - 1 - An)j 



K-5 K-3 
m=3x2^^-l,n=2 2 -1 



u{p- 



m=0,n=0 

K-3 K-5 
jf_3 -\ m=2 2 — l,n=3x2 ^ — 1 

1 - 4m) + (3 X 2— - 1 - 4n)j \ . (19) 

J m=Q,n=0 



"^opt 



The proof of Proposition [s] is given in Appendix A In principle, the optimal constellations 3^, 
and 3^op| in Proposition [s] can be attained by starting at any corner point in Z with the largest 
energy and then, for \X\ = 2, successively selecting all nearest neighbors of the previously 
already selected points along the diagonal lines and for \X\ =4, successively selecting all every 
other points of the previously already selected points along the horizontal and vertical lines. 
Fig. [T| visually demonstrates how the optimal constellations yj^^^^ and y^^^ in Proposition [s] have 
been obtained by starting at the corner point in the first quadrant with the largest energy via the 
diagonal line or the horizontal and vertical lines search over the cross QAM constellations when 
K is either even or odd. However, we should clearly point out here that the way in obtaining the 
optimal constellation 3^opt in Proposition |5] for X^pt = {1, j} is exactly the same as the one in 
partitioning a large constellation into two small sub-constellations in the trellis coded modulation 



(TCM) proposed by Urgerboeck in |79|, [80|. The major difference between these two partition 
methods is that the partition of a constellation in the TCM is based on a union operation, i.e., 
the union of the two partitioned sub-constellations is equal to the original constellation, whereas 
the partition of a constellation in the UFCP is based on multiplication. As an application of 
Proposition [5| we give the following example. 
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Example 6: = 4. By Proposition |5} we have that 
1) if \X\ = 2, then, 



(1) 



opt 

(1) 



yo^t = {3 + 3j, l+j,-l + 3j, 3-j, -3 + j,-l-j,l-3j,-3-3j}, 
2) if I A" I = 4, then, 

{l,-l,J,-j}, 
{3 + 3j,-l + 3j,3-j, 



X, 



(2) 



opt 



y opt 



IV. Collaborative Unitary UFCP Space-Time Block Codes 



In this section, we take advantage of UFCPs established in Section III to carefully design 
the constellation triple A, Bi and B2- We will show that it is the unique factorization of signals 
that results not only in the unique identification of the channel coefficients and the transmitted 
signals in the noise-free case, but in full diversity in the noise case as well. 

A. Signal Designs and Unique Identification 



In order for the code ( [3b| ) to enable the unique identification of the channel and the transmitted 
signals as well as full diversity, the three constellations A, Bi and B2 must work cooperatively. 
There may exist many cooperative agreements to make these constellations work together. 
Different collaborative ways will produce different codes, but in this paper, we require that 
A, Bi and B2 cooperate in such a way that A and Bi, A and B2 constitute two pairs of UFCPs. 
In other words, the same constellation A collaborates with both the constellations Bi and B2. 
More explanations of why such a collaborative scheme is adopted will be given after Theorem [T] 
To make the selection of these constellations more clear, we change notation A, Bi and B2 into 
the respective X, J^i and 3^2- Correspondingly, the code now takes the following form: 



U 





/ 


X 


\ 


> 


1 







X* 


: xeX^y^e yi,y2 G 3^2 > 


+ |?/i|2 + |?/2|2 




yi 


1/2 




v 


-y*2 


yl I 


> 



(20) 
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where ~ A" and 3^2 ~ constitute two UFCPs. Such a code U is called a collaborative 
unitary UFCP space-time block code. 

Now, we formally state the first main result in this paper. 

Theorem 1: Let ~ X, and y2 ^ X he two UFCPs. Then, for any given nonzero received 



signal vector without noise, i.e., r = Uh, the code U designed by pO| ) enables the unique 
identification of the channel coefficients and the transmitted signals. ■ 



Proof: Using the code structure ( [20| ), we have 

1 



Xh, (21a) 



= A 12^/ 12^ I 12 -^^' ^^^^^ 



where X and Y are defined by 



X 



X 

X* 



-y*2 y*i 

Since |3^j| > 1 for i = 1,2, by Proposition [T] we have x 7^ and as a result, the matrix X is 



invertible. Eliminating h from pT] ) yields 

r„ = YX-^r,. (22) 



X 



Notice that YX ^ = \ ^ ^ is the Alamouti codeword matrix. Therefore, equa- 

j/2 (yi\* 



tion ( [22| ) can be rewritten as 

Vy = R^.u, (23) 

where Vy = (rj^(l), r*(2))'^, u = {yi/x,y2/xY and is given by 

r,(l) -r,(2) 
r:.(2) r:(l) 

It can be verified that R^^ is unitary up to a scale. In addition, since r = (r!^, r^)^ 7^ and X is 
invertible, we have h 7^ and thus, r^^ 7^ 0, which is equivalent to the fact that R^, is invertible. 
Therefore, from ([23]) we obtain 

u = (24) 
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That is 

y, _ r:(l)r,(l)+r.(2)r;(2) 



(25a) 



X |r,(l)|2 + |r,(2)|2 ' 

y, ^ -r:(2)r,(l) + r.(l)r;(2) 

X |r.(l)|2 + |r.(2)|2 • ^ ^ 

Since X and form the two UFCPs and x E X and yi E yi for i = 1,2, once their 
quotients yi/x have been determined, then, x and y^ themselves can be uniquely determined. 



In other words, there exists a unique triple x,yi and y2 such that ([25]) is satisfied. Moreover, 
after we have determined x and yi, then, the channel vector h can be uniquely determined by 
h = a/IxP + + |?/2pX"-'-r^ = ^/\x\'^ + + \y2\'^Y^^ry. This completes the proof of 
Theorem [TJ □ 
Some observations on Theorem [T] are made as follows: 

1) Theorem [T] tells us that using the UFCP code ( [IQ] ), the channel coefficients and the 



transmitted signals can be uniquely identified by only processing the four received signals. 
Theorem 1 is very desirable in the design of noncoherent constellations, since if a signal 
design is not able to provide the unique identication of the channel and the transmitted 
signals in the noise-free case, then, the reliable estimation of the signal will not be 
guaranteed, even in high SNR. 
2) From the proof of Theorem [T| i.e., ([25]), we can observe that if either ~ A" or 3^2 ~ X, 



then, the channel coefficients and the transmitted signals can be still uniquely identified. 

In addition, ( |25a[ ) divided by ( |25b[ ) gives us 

y, _ rl{l)ry{l) + r,{2)r;{2) 

y2 -r*(2)r,(l)+r,(l)r;(2)- 
Now, it can be seen clearly that if ~ 3^2, then, yi,y2 and thus, x can be uniquely 
determined. As a result, the channel coefficients can also be uniquely determined. There- 
fore, it is worth emphasizing the fact that it is the unique factorization of constellations 



proposed in Section III that enables the channel coefficients and the transmitter signals to 
be uniquely identified. In a word, the unique identification of both the channel coefficients 
and the transmitted signals requires three constellations collaboration. 
3) The aforementioned Observations 2) tells us that there are several ways in which the three 
constellations A, Bi and B2 can collaborate so that both the channel coefficients and the 
transmitted signals are able to be uniquely identified. However, in this paper we adopt 
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yi ^ X and y2 ^ X, since the same constellation X collaborating with the other two 
constellations and 3^2 not only enables the unique identification and thus, full diversity 
(see Theorem [2]), but also provides an opportunity for the two transmitter antennas to 
accumulate their minimum Euclidean distances such that the coding gain is increased. See 
more details in Section |Vl 

B. Full Diversity 



In order to analyze the full diversity of the UFCP code pO] ), let us first consider a general 
space-time block coded noncoherent MIMO system with M transmitter antennas, N receiver 
antennas and flat fading channels as follows: 

T = SH + H, (26) 

where T denotes a T x N received signal matrix, H denotes an M x channel matrix, 
S is a T X M codeword matrix and H denotes a T x noise matrix. We assume that the 
elements ^tn of H are samples of independent circularly- symmetric zero-mean complex Gaussian 
random variables with variance a^. Under these assumptions, the probability density function 
of the received signal matrix T conditioned on the transmitted signal matrix S is the Gaussian 
distribution, i.e., 

1 / Tr(T^(SS^ + (T2l)-iT)' 



X exp 



7r™det(SS^ + a2I) ^ \ 

and thus, its likelihood is given by 

Tr(T^(SS^ + (T2l)-iT) 



cr2 



- In det(SS^ + a^I) - TiV In TT. 



Then, the maximum likelihood (ML) receiver for the noncoherent MIMO system is equivalent 
to solving the following optimization problem: S = argminsjA^ — Tr(T^©sf)}» where 
= (a^I + S^S)"^ and = Indet (a^I + S^S). To avoid estimating the variance 
of noise, consider the conditional probability density function of the received signal matrix T 

given the channel matrix H and the transmitted signal matrix S, i.e., 

12 ~ 
If 



1 / ||T-SH|i2 



X exp 



and thus, its likelihood is given by 



1^ ^^Ik _TArin7r-2TiVlna. 



Ct2 
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The generalized likelihood ratio test (GLRT) receiver for the joint estimation of H and S is to 
maximize the likelihood, which is essentially equivalent to solving the following nonlinear least 



square error optimization problem |7TJ, [|8T|, [82| 



{H,S} = argmin||T - SH||^. (27) 

H.S 

Its solution can be obtained by first estimating the transmitted signal matrix S as 

S = arg max Tr (t^S (S^S) S^t) , (28) 

and then, estimating the channel matrix H as H = (S^S) ^S^T. Particularly for any unitary 
code, the ML receiver and the GLRT receiver for the optimal estimation of the transmitted signal 
matrix are equivalent, i.e., S = argmaxs Tr (T^SS^T). In addition, Brehler and Varanasi |7l| 
analyzed the asymptotic performance of the GLRT detector for the noncoherent MIMO system 
and proved the following lemma. 

Lemma 1: Let a 2M x 2M matrix R^^ be defined as H^g = (S,S)^(S,S). If each matrix 
Uss has full rank for all pairs of distinct codewords S and S, then, the resulting space-time block 
code provides full diversity for the GLRT receiver, and moreover, the pair-wise error probability 
Pglrt(S — 7- S) of transmitting S and deciding in favor of S 7^ S has the following asymptotic 
formula: 

f 2MN - 1 \ . ^ . 

det^(S^S) 

Pglrt(S ^ S) = ^ r^vT^T^ X SNR-*^^ + o(SNR-^^). 

det [rlssj 

■ 

Lemma [l] tells us that the full rank of the matrices R^s for all the distinct codewords S and S 
assures full diversity. It is not difficult to prove that a necessary condition for R^^ to have full 
rank is T > 2M. Hence, in this paper, we consider the case of the shortest coherent time slots 
when it is possible for the UFCP code design to enables unique identification of both the channel 
coefficients and the transmitted signals as well as full diversity. Now, we are in a position to 
state the second main result in this paper. 

Theorem 2: Let ~ A", and y2 ^ X constitute two UFCPs. Then, the code U designed 



by (20) enables full diversity for the noncoherent ML receiver. Furthermore, the pair- wise error 
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probability Pml(U — > U) of transmitting U and deciding in favor of U 7^ U has the following 
asymptotic formula: 

3 



Pml(U ^ U) 



X 



+ o(SNR-2), 



where R„^i = (U, U) " (U, U) . ■ 

Proof: We first note that since the UFCP code (|20j) is unitary, we have det(U^U) = 
det(U^U) = 1. By Lemma [l| we only need to prove that Huii is invertible for any pair of 
distinct U and U. Since (U, U) is a square matrix and Kuu = (U, U)^(U, U), proving that 
the matrix Huu is invertible is equivalent to proving that the matrix (U, U) is invertible. Notice 
that 

/ 



(u,u) 



a/|x|2 + ||/i|2 + |y2|2 



X 





X 


\ 





X* 





X* 


yi 


1/2 


yi 


h 


-y2 


yi 


-2/2 


yl ) 



(29) 



By some algebraic manipulations, we can obtain that the determinant of (U, U) is given by 



det(U,U) = (It 



yi 



+ 



y2 



y2 



+ I2/1P + 1^2 n • Since X and 3^, for i = 1,2 
forms two UFCPs for any U 7^ U, i.e., (x,?/i,?/2) 7^ {x,yi,y2), we have det (U, U) 7^ 0. In 
addition, the asymptotic formula can be immediately obtained by utilizing Lemma [T] again with 
M = 2, N = 1 and S = U. This completes the proof of Theorem [2} □ 
From the proof of Theorem |2] we can observe that it is the unique factorization of our designed 
constellations that enables the matrix Hyy to have full rank. In addition, we can also observe 
that the condition in Comment 2) on Theorem [T] is still a sufficient condition for the matrix 
Hyy to have full rank. Therefore, noncoherent full diversity also requires the three constellations 
collaboration. 



V. Optimal Designs of Unitary UFCP Stace-Time Block Codes 

Our main task in this section is to efficiently and effectively optimize the coding gain for 
the unitary UFCP space-time block codes generated from the energy-efficient cross QAM 
constellations. 
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A. Problem Formulation 

Theorems [T] and [2] together tell us that the unitary UFCP code designed by pO] ) enables 
the unique identification of the channel coefficients and the transmitted signals as well as full 
diversity for the noncoherent ML receiver. Therefore, the code design partially gives a solution 
to Problem [T] In order to further optimize its error performance, we can see from the asymptotic 
formula of the pair-wise error probability in Theorem [2] that when SNR is large, the error 
performance is dominated by the term det ((U, U)^(U, U)). Hence, following the way similar 
to coherent MIMO communications [4|, we define the coding gain for the unitary code as 



mm A/det ((U,U)^(U,U)). 



(30) 



Theoretically speaking, we should maximize the coding gain G'(A:', 3^i, 3^2) directly among all 
two UFCPs yi ^ X and 3^2 ~ X. However, this optimization problem, in general, is too 
difficult to be solved, since the optimal design of constellations is generally extremely difficult 
to be reformulated into a tractable optimization problem, even for an additive white Gaussian 



noise channel [78 1, [83|-[88|. To make the problem tractable, in this paper, we restrict ourselves 
to using the energy-efficient cross QAM constellations to generate the two UFCPs yi ^ X and 
3^2 ~ with X C { — 1, —j, 1, j}. Specifically, let Zi and Z2 be two given 2P-ary and 2'^-ary 
cross QAM constellations, respectively. Then, the three constellations X, 3^i and 3^2 are selected 
in such a way that Zi = ^ and Z2 = In addition, it is not difficult to verify that if Z- — — 
and ^2 = 37, then, aZi 



ay 



- and aZo 



X 

— for any positive a. Therefore, a family of UFCP 



codes resulting from the cross QAM constellations and an energy scale a is characterized by 



W„(A',3^i,3^2) = <U„ 



1 



X 


o\ 


> 







X* 


: a; G A*,?/! G 3^i,?/2 e 3^2 


> .(31) 


«2/i 


"2/2 






-«Z/2 


«2/7 


> 





By employing the code structure (|3Tj) and performing some algebraic manipulations, the 



expression ( [30| ) can be further simplified into 

G„(A',3^i,3^2 



y\ 12 



+ 



yi 12 



mm 



I — P 



+ 



1^2 12 



min Gc,{x,y;x,-y\X,yi,y2) 

{x,yi,y2)^ix,yi,y2),x,xGX,yi,yi<^yi,y2,y2ey2 



))(l + a2(|||2+|||2)) 

(32) 
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where notation Ga{x,y; x,y\X ,yi,y2) is defined as 



X ' ' X X 



GJx, y; X, y\X, 3^i, 3^2) = ^ ^ , (33) 

aK ,J| ,^1,^2; (^^^2(|^|2 + |^|2))(l + cfc2(|^|2+ |||2))' 

which is called a coding gain function. Our design problem is now formally stated as follows: 
Problem 2: Let |Wq,(A:', 3^i, 3^2)! = 2*^ (r > 4) be fixed. Find an energy scale a, three 
nonnegative integers 5, p and g satisfying a total transmission bits constraint: p + q — 5 = r, and 
the unique factorizations of a pair of the 2P-ary and 2'^-ary cross QAM constellations Zi and 
Z2, i.e., Zi = yi/X and Z2 = y^l ?C with X C 1, j} and \X\ = 2\ such that the 

coding gain G'q(A:', 3^i, 3^2) is maximized, i.e., 

|5,5,p, g, A',3^1,3^2} = arg max max max GafA", 3^i, 3^2)- 

^ p+q-S=rZi=yi/X,Z2=y2/X a 



B. The Solution to Problem |2] 

In order to solve each individual optimization problem in Problem |2} we first introduce some 
notation for discussion simplicity. Recall that Q denotes the modified 2^ -ary cross QAM defined 
in Definition |4j Let P denote one of its comer points with the largest energy E, and let Pi and 
P2 denote the two nearest neighbors of this corner point, with the respective energies and Et, 
where E^ > Et. For given positive integers u, v and w satisfying u + v = w with u > v, let Qi 
and Q2 denote the respective modified 2"-ary and 2''-ary cross QAM constellations; denotes 
one of the corner point in Qi with the largest energy Ei; Its two nearest neighbors are denoted 
by Zii and Zi2, respectively, with energies being En and Ei2, where En > Ei2. Specifically, for 
a given positive integer w > A, integers u and v are defined as follows: 

u = v = ^ if w is even, 



(34) 

- if 7/) iQ nrlrl 

2 ' 2 



u = i^,v = ^ if wis odd 



Correspondingly, all notations Qi.En and £',2 are defined in the same way as Qi, En and Ei2. 
As will be seen shortly in the following discussions on the solution to Problem [2} the integers 
u, V and w are to be replaced by p, q and r — 6, respectively, and thus, Qi is to be regarded as 
Zi. Hence, the corresponding positive integers p and q and energy notation eI^\ e\^ and e'i2 
are used directly without the need of redefinitions. Some properties regrading these energies are 
collected as the following Lemma. 
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Lemma 2: For the modified 2^''-ary cross QAM Q, the following statements are true: 

1) li K = 3, then, each comer point has only one nearest neighbor, = 10 and = 2. 

2) If K is even, then, E = 2(2f - l)^ and E^ = Et = (2f - l)^ + (2^ - 3)^ 

3) If K is odd and greater than 3, then, E = (2^ - 1)^ + (3 x 2^ - 1)^, 

E, = (2^ - 3)2 + (3x2^- 1)2 and E^ = (2^ - 1)^ + (3x2^- 3)^. ■ 
Lemma [2] can be verified directly by the calculation and thus, its proof is omitted. Now, applying 
Lemma |2] to a pair of QAM constellations Qi and Q2 yields Lemma |3] 

Lemma 3: Let Qi and Q2 denote the modified 2"-ary and 2^-ary QAM constellations, 
respectively, with u > v > 2. Then, 

E2 + En < E1 + E2U (35a) 
E2 + E12 < Ei + E22. (35b) 

■ 

The proof of Lemma |3] is provided in Appdendix |Bj 

Lemma 4: Let Qi and Q2 denote the modified 2"-ary and 2^-ary QAM constellations, 
respectively, with u > v > 2. Then, the following four inequalities hold: 

+ > E1 + E2, (36a) 



V^i + E21 + + E22 > V ^1 + ^21 + V ^1 + ^22, (36b) 



V^i + E2 + ^Ei + E21 > JEi + E2 + yEi + E2i, (36c) 



V^i + £^21 + VE2 + Eu > \/ El + E21 + \/ E2 + Eu. (36d) 

■ 

The proof of Lemma |4] is postponed to Appendix |C} We also need the following lemma, whose 
proof is given in Appendix |D} 

Lemma 5: Let Qi and Q2 represent the modified 2"-ary and 2''-ary QAM constellations with 



u and V defined in ( [34| ). Then, 

El + E22 < E2 + Ell, when w is even, (37a) 
E1 + E22 > E2 + E11, when w is odd. (37b) 



September 2L 2011 



Submitted to IEEE Trans, on Information Theory 



XIA, ZHANG AND DUMITRESCU, ENERGY-EFFICIENT FULL DIVERSITY COLLABORATIVE UNITARY CODE DESIGNS 



24 



Now, it is time to start solving Problem [2j Let us consider the following three cases of all the 
possible values of 6. 

1) 5 = 0: In this case, the constellation X includes only one element, i.e., X = {x}, where 
X = ±1, ±j. Since r = Uh + ^ and matrix diag(x*l2, 12) is unitary, we have diag(a;*l2, 12)1" = 
diag(x*l2,l2)Uh + diag(a;*l2,l2)^, with diag(a;*l2, 12)^ having the same statistical property as 
$, and 

/ 1 \ 

1 



diag(x*l2,l2)Uc 



(38) 



\ -ay 2 "l/i / 

Therefore, we only need to consider the case where x = 1. In this situation, ^1 = 3^1,^2 = 3^2 
and Ga{x,y;x,y\X,yi,y2) is reduced to 



G'a(l,y; l,y|^,3^i,3^2) = 7^ — 2/, ,2 , 1 i2\\f^ , — 2f\- 12 , 1- 12^^ - ^^^^ 

(1 + a^{\yi\^ + \y2n){l + a^{\yi\^ + 1^2^)) 

Let us first consider a special case where \Z2\ =4. In this case, \y2\ = \y2\ = 2 and hence, 



It is very interesting to observe that for any fixed a, X and 3^j, the objection function 
GQ,(a;, y; a;, y)|'^, 3^1, 3^2) with respect to the variables x,x,'y and y is minimized when its 
numerator achieves the minimum and simultaneously, its denominator achieves the maximum, 
both optimums being achieved when y2 is the nearest neighbor of y2 and yi = yi = Zi is the 
comer point with the largest energy Ei in Zi. Therefore, the minimum in this case is given by 

GalA", 3^1, 3^2) = min G„(l, y; 1, y\X , 3^i, 3^2) = 4a/(l + a(2 + E^)f. (41) 

Now, we consider a general case where \Z2\ > 8. In this case, notice that 

G'q(A:',3^i,3^2) = niin{ min y; x, yjA", 3^i, 3^2), 

(lmM?/2|)={|yi|,|j/2|),either|yi|=v^or|y2|=v^ 

min G'Q(a;,y;x,y|A',3^i,3^2), 

(ll/lMl/2|) = (|yi|,|i/2|),|s/l|>V^,||/2|>V2 

min G'a(x,y;x,y|A:',3^i,3^2)}- (42) 



Realizing the following facts is key to obtaining the minimum of the objection function ( [39| ) 
with respect to variables x, x, y and y: 



September 2L 2011 



Submitted to IEEE Trans, on Information Theory 



XIA, ZHANG AND DUMITRESCU, ENERGY-EFFICIENT FULL DIVERSITY COLLABORATIVE UNITARY CODE DESIGNS 



25 



(a) Under the conditions that y 7^ y, \y2\) = {\yi\, Iml) and either \yi\ = a/2 or jyal 
\/2, the numerator of the objection function (|39l) is lower-bounded by Aa"^, i.e., 



aWvi - yi\^ + \y2 - ml"^) >^oi^ 



(43a) 



where the equality holds when either ?/2 is the nearest neighbor of y2, i.e., |t/2 — ^2! = 2, 
and yi = yi or yi is the nearest neighbor of yi and ?/2 = ^2- The denominator of the 



objection function p9) is upper-bounded by 

(1 + a\\y^\' + |y2p))(l + a\\m\' + Iml')) < (1 + + 2))^ 



(43b) 



where the equality holds when both yi and yi are the comer points in Zi with each having 
the largest energy Ei. Therefore, both equalities in ( |43a| ) and ( |43b[ ) hold simultaneously 
when both ?/i and yi are the comer points in Zi with each having the largest energy Ei 
and y2 is the nearest neighbor of ^2- 
(b) Under the conditions that |y2|) = {\yi\, Ihl), Ivil > v^, 1^/21 > and y 7^ y, the 



numerator of the objection function (39) is lower-bounded by 



Qi^(i?/i -yip + 12/2 -y2n > 



if \yi\ > En and \y2\ > E21, 
if either \yi\ < En or I1/2I < -f^: 



(44a) 



The denominator of the objection function ([39]) is upper-bounded by 

(l + a\\y,\' + \y2mi + a'{\m\' + \y2\')) 
[1 + a\Ei + E2)f if \yi\ > En and 1^2! > ^21, 
;i + a\Ei + E2i)f if either \yi\ < En or |?/2| < ^21- 



< 



(44b) 



(c) If \Zi\ > IZ2I and {\yi\, I1/2I) 7^ (|yi|, 1^21), then, Ei > E2. Note that the numerator of the 



objection function ( [39| ) is lower-bounded by 



aWyi ~ yi\^ + \y2 - ml^) > 



(45a) 



where the equality holds when either y2 is the nearest neighbor of y2 and yi = yi or yi 
is the nearest neighbor of yi and ?/2 = ^2- In addition, the denominator of the objection 



function ( [39| ) is upper-bounded by 

(1 + + |y2r))(l + a'dyiT + 1^2!')) < (1 + a\Ei + E2)){1 + ^^(Ei + E21)), 



(45b) 
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where the equality holds when both yi and yi are the corner points in Zi with each 
having the largest energy Ei, while one of y2 and y2 is the corner point in Z2 with the 
largest energy E2, and the other is of the second largest energy £'21- In other words, the 
upper bound (1 + a'^{Ei + E2)){1 + oP'{Ei + E21)) is the second largest maximum of the 
denominator of the objection function (|39l). Very interestingly, under this condition, both 



equalities in ( |45a[ ) and ( |45b[ ) are able to be achieved at the same time when yi = yi is the 
comer point having the largest energy Ei in Zi, y2 is the nearest neighbor of y2 and of 
the largest energy E2 and y2 has the second largest energy £'21 in Z2. 
The above three observations reveal that 



min G'a(l,y; l,y|A',3^i,3^2) 

i\yi\,\y2\)={\m\,\y2\),\vi\=V2or\y2\=V2 

4a2 



mm 

i\yi\,\y2\)={\yil\y2\),\yi\>V2,\y2\>V2 



{l + a^{E, + 2)y' 
Go.{l,y-l,y\X,y,,y2) 



(46a) 



> 



[l + a^{E, + E2i)r 



(46b) 



min Ga{l,y;l,y\X,yi,y2) 

{\yi\,\y2m\yi\,\y2\) 

_ 4^2 

" (1 + + E2)){1 + + E21)) ■ 

When 1^2 1 > 4, Ei > E21 > 2 and as a result, comparing ( |46a| ) with ( |46c| ) gives us 

4^2 4a2 

{l + a^{Ei + 2)y ^ (1 + a2(Ei + E2))(l + ^2(^1 + E21)) 

and thus, we have 

min Ga{x,y;x,y\X,yi,y2) 

{\yi\,\y2\) = {\yi\,\y2\),\yi\=V2or\y2\=V2 

> min Ga{x,y;x,y\X,yi,y2). 

i\yi\,\y2m\yil\y2\) 



(46c) 



(47a) 



In addition, since E2 > E21, we can obtain from (46b) and (46c) that 



min Gc,{x,y]x,y\X,yi,y2) 
{\yi\,\y2\)=(\yil\y2\),\yi\>V2,\y2\>V2 

{\yi\,\y2\myi\,\y2\) '^^'^2; 



(47b) 
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Combining (|42]) with (|47]) yields 



Ga{X,yi,y2) = —— jT-= „ XX , JT-^ , • (48) 

(1 + ^^(Ei + -£/2))(l + a^{Ei + E21)) 



Since if IZ2I = 4, i?2 = -£^21, the conclusion d48j) includes d4Tb as a special case. Therefore, no 



matter whether or not IZ2I > 4, we always have ( [48] ). Notice that ( [48] ) can be rewritten as 

^,(^,3^1,3^2) = 



4 



(a-i + + E2))(a-i + a{Ei + E21)) 

4 

a-2 + a2(Ei + E2)(Ei + E21) + {2Ei + E2 + E21] 



(49a) 



< — ^ (49b) 

2v/(^i + ^2)(^i + ^2i) + (2^1 + ^2 + ^21) 

4 



(a/ -El + E2 + v^-Ei + -£^21) 



2 ' 



where we have used the arithmetic mean and the geometric mean inequality from ( |49a| ) to ( |49b[ ): 

a + b > 2\/ab, and the equality holds when a = 1/ \J {E\ + E2){E\ + -^21)- Thus, we have 

4 

maxG'Q(A',3^i,3^2) = , ,2 

(V^i + E2 + V^i + E21) 

4 

< , , (50) 



where we have used Lemma |4] with w = r^u = p,v = q and Qj = and the inequality in ( |50l ) 
holds when p = p and q = q. All the above discussions can be summarized as Property [T] 
Property 1: When p > q and 5 = 0, the optimal solution to Problem |2] is given as follows: 

(1) If r is even, then, p = q = r /2. 

(2) If r is odd, then, p = (r + l)/2, g = (r - 1) /2. 

Once the optimal p and g have been determined, then, X^^^ = {1}, y[^^ is the 2^-ary QAM 
constellation Zi, 3^2°^ is the 2''^-ary QAM constellation Z2 and the optimal energy scale a is 



determined by 

5= , ^ (51) 

^(E{n^f)(^;n^(?) 

Moreover, the optimal coding gain is given by 

yi'\ yi'^) = \ (52) 
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In fact, Property [T] gives us an optimal design of the unitary training Alamouti code based on 
the cross QAM constellations. 

2) 5 = 1: In this case, there are in total six possibilities in choosing the constellation X, i.e., 

X = {1, -1}, X = {1, j}, X = {1, -j}, X = {-1, j}, X = {-1, -j}, X = {j, -j}. In order to 
attain the optimal solution, we take the following two steps. 

Step 1: Determine G^{XW ^yl'lyi''^) = min(,,yT)^(£,yT) y; x, y| 3^J^\ 3^^^)), 

where constellations X^^^ and y^^^ are derived from Proposition [s] with Z = Zi, i.e., X^^^ = 
{l,j} and = 3^i,opt. To do that, we first represent Gc,{X^^\y[^\y^^^^) as 

Ga{X^'\yi'\yi'^)=mm\ min G^{x,y;x,y\X^'\y['\yi'^), 

min G^ix,y;x,y\X^'\yi'\yi'^)], (53) 

which leads us to individually considering the following two optimization problems, 
(a) a; = X. In this case, the objection function in ([32]) is simplified into 

r vr vi;p« v« v^'h - ^Hlyi - ml^ + \y2 - ml^) .... 

(1 + a^{\yi\^ + |i/2r))(i + a^imr + mn) 

since |x| = = 1. Following the discussion very similar to the case when 5 = 0, we can 
obtain 

min G^{x,y;x,y\X^'\yi'\yi'^) 
(^,y j'Fi^^y ),x=x 



(l+a2(Si+£;2))2 
8o" 



(l+a^{Ei+E2)){l+a2{Ei+E2i)) 



ifg = 5org = 2org = 4 and p = 5, 
if g = 3, 



if q is even and greater than 2 and p ^ 5 
I ,1 , 2117 , TP I I F ^^ if Q is an odd integer exceeding 5. 

{l+a^{Ei+E2i)j{l+a^(Ei+E22)) ^ to to 



(55) 



{l+aHEi+E2i))^ 



(b) X ^ X. Under this condition, we further split the feasible domain of the optimization 
problem: min(j. yT)-^(£ yT) ^..^^ Ga{x, y; x, y\X^'^\ 3^j^\ 3^2^^), into four disjoint sub-domians 
as follows: 

1^11 = {(a:,y^x,y^) : x ^ x, {\y,\, \y,\) = (|y2|, 1^21) = (^2,^2)}, 

= {(x,y^,x,y^) -.x^x, {\y,\, \y,\) = iE,,E,), (ly^l Iml) 7^ (^2,^2)}, 
1^13 = {(x,y^x,y^) : X ^ X, ^ (^1,^1), (I2/2I, Ihl) = (^2,^2)}, 
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2^14 = {(x,y^,x,y^) : x^x,{\y^l\y^\) ^ (Ei, Ei), (l^/al, 1^21) ^ (^2,^2)}. 
Therefore, we have 



mm 



(56) 



Now, let us first consider each inner minimization problem, 
a) When (x,y^,x,y^) G Dn U I?i2 and |Zi| > 8, |^ - f | > ^20 and as a result, the 



numerator of Ga{x,y\x,y\X''^\y']^\y':^^) is lower bounded by 



> 2W. 



yi 


_ m 


2 

+ 


m 


_ m 


X 


X 




X 


X 



(57) 



Under the same condition, the denominator of Ga{x, y; x, y|A:'(^\ 3^i^\ 3^2^"*) upper 
bounded by 



|2 ^ 1^212 



1 



1^1 |2 ^ 1^212 



X 



X 



Combining ( [57| ) with ( |58] ) results in 

G«(x,y;x,y|A'«,3^«,3^^^))> 



mm 

(x,y'J^,£,yT)eDiiUX'i2 



(58) 



(59) 



[l + a^{Ei + E2W 

b) When (x, y^,x,y^) G ©14, we first notice that the numerator of the objective is 
lower bounded by 

2 



yi 


_ m_ 


2 
+ 


y2 


_ m 


X 


X 




X 


X 



> 8a' 



(60) 



since I ^ - ^ I > 4 and 

\ X X \ — 

is upper bounded by 



y2 y2 



> 4. In addition, the denominator of the objective 



|2 _^ 1^2 |2 



X 



X 



1 + 



yil2 _^ 1^2 |2 
X ' ' X ' 



< (1 + a\E, + E2i)){l + a\E,^ + E^)). 



(61) 



Now, let us argue that both equalities can be achieved simultaneously. Since x, x G 
X = {1, j}, without loss of generality, we can always assume that x = 1 and x = j. 
Recall the definition and properties of the Group-x which we have discussed in 



Section III If we use 1 and Zij to denote two groups generated by J^j ~ A* for 
i = 1, 2, then, we have Zi i U Zij = Zi and Zi^i fl Zij = $. Also, recall that notation 
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Zi denotes one of the corner points in Zi with the largest energy Ei, and notation Zix 
and Zi2 denote its two nearest neighbors in Zi, i.e., \Zi — Zn] = \Zi — Zi2\ = 2, with 
the respective energies En and Ei2. Then, Zn E Zij for i = 1,2, since i = 
and (imin(3^j) > 2. Hence, if we set i/i = Zi and fji = jZn, then, both the equalities 
in (|60l) and dMl) hold at the same time and thus. 



min G^(x,y;x,y\X^'\yi'\y^2'^] 



(62) 



(1 + a^iE, + E2i))(l + a^iEn + E^)) ' 
c) When (x,y-^,x, y^) G I'ls, we need to consider two possibilities: 1 = 4 and 
1 > 8. If 1^2 1 > 8, then, following the discussion similar to Situation a), we 
obtain 

min G.ix,y;x,y\X^'\yi'\y^'^) > ^ (63) 

If IZ2I = 4, then, \y2\ = \y2\ = 2 and li/2/a; — ^2/2^! = 2. Following the same 
discussion as Situation b) and choosing yi = Zi,yi = jZu, we have 

min Go.{x,y;x,y\X^'\yi'\yi'^) 

{x,y^ ,x,y^ )eI5i3 

8«2 

(64) 



{l + a^{E^ + 2)){l + a^{En + 2))' 
If we make the convention that E2 = E21 = 2 when Z2 is the 4-QAM constellation, then. 



equation ( [62| ) includes ( [64| ) as a special case. To compare ( [59| ) with ( [62| ), we need to prove 
that the following inequality 

3(1 + a\Ei + E2i)){l + a^{En + ^2)) > (1 + + ^2))' (65) 
is true for any positive a. To do that, we first establish an energy inequality: 

3^1^11 + 3EiE2 + 3^11^21 + 3^2^21 > (^1 + E2)^. (66) 

To show this, let us discuss the following four possibilities, 
a) \Zi\ = IZ2I = 8. In this special case, it can be verified by calculation that 3-E'i-E'ii + 



3^1 ^2 + 3^11^21 + 3^2^21 = 432 and {E^ + E2Y = 400. Thus, the inequality §E\ 
is true in this particular situation. 
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b) IZ2I = 4. Since Si^n +A>Ei, we have 



C) IZil > 



SEiEii + 3E1E2 + 3i?iii?2i + SE2E21 
= E^{3En + 4) + 2(3En + 4 - Ei) + 4Ei + 4 
> ^2 + 2^1^2 + 4 = (El + ^2)=^. (67a) 
8, In this case, note that 2Eii > Ei,Ei > En > E2. Hence, we attain 



SEiEii + 3E1E2 + 3i?ii-E'2i + SE2E21 
> 1E\E-ii + 2E1E2 + EiEii 
= E\ + 2E1E2 + -^2 ~ (-^1 + -^2)^- 
d) iZil > 12:21 > 8. Then, 2Eu > Ei and 2^21 > ^2 and as a result. 



(67b) 



SEiEii + 3E1E2 + 3-Eii-E'2i + SE2E21 

> 2E1E11 + 2E1E2 + 2E2E21 

> El + 2E1E2 + ^2' = (^1 + ^2)'. (67c) 
Combining ( |67a| ), ( |67b[ ) with ( |67c[ ) gives us 

3(1 + a2(Ei + E2i))(l + a'(^ii + ^2)) 

= 3 + 3(Ei + ^2 + Ell + E2i)a^ + (3EiEii + 3E1E2 + 3E11E21 + 3E2E2i))a^ 
> 1 + 2(Ei + E2)a^ + (El + E2)2a^ 

= (l + a(Ei + E2))l (68) 
Hence, the inequality ( [65] ) is true. Now, comparing ([59]) with ([62]) and using ( [65] ) lead to 

min(^,yT^£^yT)g2?i4 ^^(x, y; x, y| A'(i), yf \ y^^^) 
. 2(l + a^(Ei + E2i))(l + a^(Eii + E2) 

- (l + a^(Ei + E2))^ ^''^ 

if 1^2! > 8, which is equivalent to the fact that 



G^ix,r,x,y\x^'\yi'\y^2'^) 

> min G^ix,y;x,y\X^'\y['\yi'^). 



min 

(x,y'J^,£,yT)eI'iiUX'i2 



(70) 
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Using the same argument, we can derive that 



min G^{x,y-x,y\X^^\y\'\y^^^) 

(x,y^,x,y^)el'i4 

< min G„(x,y;x,y|A'(i),:yf\3^(')). 

(x,y^,a:,yT)eX'i3 



(71) 



Combining (56) with (70) and (71 ) altogether tells us that 



min G^{x,y;x,y\X^^\y\'\y^^^) 

{x,y^ ,x,y^ ),x^x 



~ (1 + a^Ei + ^2i))(l + a2(En + E2)) ' 
Substituting (|55]) and (|72]) into ([53]) yields 



(72) 



(iWCBi+Ba))^ 
8a2 



if p > g = 5 or p = 5, g = 4, 
if p > g = 3, 



(l+a'2{Ei+E2)){l+a-^{Ei+E2i)) 
{l+a-^{Ei+E2i))(l+cfi{E^+E22)) P > 7^ 3, 4, 5 

■v2i'R, _l-E^,^V1-l-^2|'E,,4_E,^^ if p = g 7^ 3, 5. 



(73) 



(l+a'^{Ei+E2i)){l+a'^{Eii+E2)) 

Step 2: Establish the achievable upper bound of the coding gain for any fixed UFCP, i.e., 



G^{x,yi,y2)<G^{x^^\y\'\y!i^) 



(74) 



for any positive a. As we have mentioned before, when 5=1, there are totally six candidates 
regarding the constellation X, i.e., X = {1,-1}, A* = {l,j},X = {l,—j},X = {—l,j},X = 
{—l,—j},X = Here we only consider the case where X = {1,-1}, since the 

discussion for the other cases are very similar. Now, we examine the following possibilities: 

Case 1: g = 3. In this case, since Z2 includes the 4-QAM constellation as a subset, there are 
two points 1/20 and ^20 in 3^2 such that |?/2o - ^20! = 2 and |?/2o| = |^2o| = 2 and thus, 

GaiX,yi,y2) = min ^^(a;, y; x, y jA", 3^i, 3^2) 

< Ga{{l,Zi,y2o;l,Zi,y2o)\X,yi,y2) 

4a2 

(l + a2(E^ + 2))(l + a2(^^ + 2))- 

Case 2: g = 5. There are totally 8 corner points in Z2 with the largest energy and 3^2 

includes four of them. If two of the four points, say 1/20 and y2o, are in the same quadrant, i.e., 

\y2o -y2o\ = 2^2, then, we have 

Ga{X,yi,y2) = min Ga{x,y;x,y\X,yi,y2) 
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< Ga{l,Zi,y2o;l,Zi,y2o\X,yi,y2) 
8«2 



Case 3: p > q and q ^ 3,5. 

(a) Z21 e ^2- Then, we have 

Ga(X,y,,y2) - min Gc,{x,y;x,y\X,yi,y2) 

< Ga{{i,z^,Z2;i,z,,Z2i)\x,yi,y2) 

(1 + a^{Ei + E2i))(l + a^{Ei + ^2)) ' 

(b) Z22 e ^2- Similarly, we can obtain 

Ga{X,y^,y2) = min Ga(x,y;x,y|A',3^i,:^2) 

(a;,y^)^(£,y^) 

< Ga{i,z,,Z2;i,z,,Z22)\x,yi,y2) 

(1 + «2(^^ + ^22))(1 + + ^2)) ■ 

(c) Z21, Z22 ^ 3^2- Then, — Z21, — Z22 € 3^2 and as a consequence, we arrive at 

Ga{X,yi,y2) = min Gaix,y;x,y\X,y,,y2) 

< ^^((i, Zi, -Z21; 1, Zi, -Z22) I A", 3^1, y2) 

8a^ 



(75b) 



(75c) 



(75d) 



(75e) 



(1 + a^{Ei + ^2i))(l + + ^22)) ■ 

Case 4: p = q and q ^ 3,5. We consider the following possibilities: 

(a) If either Zu e yi or Z21 G 3^2, then, following Case 3-(a), we can have 

Gq,(A',3^i,3^2) < TT— — ^7^^ — , NW-, , — ^7^^ — , NX ) ('75f) 
(1 + a-'^Ei + -£/2i))(1 + a^{Ei + E2)) 

where Ei — E2 and E^n = £^21- 

(b) If either Z12 e 3^i or Z22 £ 3^2> then, similar to Case 3-(b), we can obtain 

G^{A:,y^,y2) < ^^ + ^2^E, + E22)){l + a-{E, + E2)y ^^^^^ 
where Ei — E2 and £'12 = -E'22- 

(c) Zu,Zi2 i 3^1 and ^21,^22 i 3^2- Then, -Zn,-Zi2 G 3^1,-^21,-^22 e 3^2 and as a 
result, we attain 

G'a(A',3^i,3^2) = min G^{x,y;x,y\X,y^,y2) 
(a;,y^)^(£,y^) 
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< ^,((1, Zi, -Z21; -1, -Zn, ^2)1^, 3^1, 3^2) 

„2 



(1 + + E2i))(l + + Eu)) ■ 



(75h) 



Comparing ( [75] ) with ( [53| ) gives ( [74| ) as required. Now, using the geometrical and arithmetical 
mean inequality and then, applying Lemma |4] to (|74]) and ( |73| ) result in 

where 



(76) 









if p > g = 5 or p = 5, g = 4, 








if p > g = 3, 














if p > g,g 7^ 3,4,5, 














if p = g 7^ 3, 5. 






,Ei^^y 



(77) 



All the above discussions can be concluded as Property |2] 

Property 2: When 5 = 1, one of the optimal solution to Problem [2] is that X^^^ = {l,j} and 



that p, q and y^^^ are determined as follows: 



(1) If r is even, then, p = {r + 2)/2,q = r/2, ^ = y 



i,opt 



and 



a 



^iE[^hEg^mE[^^+Eg^) 



if r = 10 or r 

if r = 6, 

if r ^ 6,8,10. 



(78) 



Furthermore, the optimal coding gain is given by 



Ga{x,yuy2] 



E['^+Ei'^ 



5(l)_L?(l)a_./g(l)_L?{l) 



{JE\^>+E 



E\''+Ei,{')2 



if r = 10 or r 
if r = 6, 

if r ^6,8, 10. 



(79) 



(2) If r is odd, then, p = g = (r + l)/2, ^ = y, 



i,opt 



and 



(1) 



a = < 



if r = 9, 
if r = 5, 
if r ^ 5,9. 



(80) 
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Furthermore, the optimal coding gain is given by 



if r - 9, 



if r ^ 5,9. 



Actually, Property [2] also suggests us that for 5 = 1, the optimal UFCP designed by Proposition [5] 
is still optimal in the sense of maximizing the coding gain. 

3)5 = 2: In this case, X = {1, — 1, j, — j}. We examine the following possibihties: 

(a) If either Zu G or Z21 E 3^2^ then, we can have either 

Ga{X,yi,y2) = min Ga{x,y;x,y\X,yi,y2) 
(a;,y^)^(£,y^) 



or 



< G^{l,Z^,Z2■,l,Z^^,Z2\X,y^,y2) 



(1 + a^{Ei + E2)){1 + a^{En + E2)) ' 



Go.{X,y^,y2) < G,(l,Zi,Z2i;l,Zi,Z2|A',3^i,3^2) 

4«2 



(82a) 



(1 + a^E^ + E2i))(l + a2(Ei + E2)) ■ 
(b) If either Z12 G or Z22 G 3^2, then, similar to situation (a), we can obtain 

4a2 



(82b) 



G'a(A',3^i,3^2) < max{— — — — — 

(1 + a^[Ei + E2))[l + a^[Ei2 + E2)) 

4^2 ^ 



^1 + a^{Ei + E22))(l + + E2)) ■ 

4a2 



(82c) 



(1 + a2(^^ + E2)){1 + a2(Ei2 + ^2)) ' 
since E12 + E2 < Ei + E22 according to Lemma [3| 
(c) If Zn, Z,2 ^ yi and Z21, Z22 ^ 3^2, then, Zn, ^12 G 2i,-i U Zi,, U and Z21, Z22 G 

^2 -1 U Z2J U ^2 - J- In this case, one of the following two possibilities must occur: 
a) If either Zu, Z\2 G ^i,zi or Z21, Z22 G 22,X2 5 then, we have either 

G'„(A',3^i,3^2) = min y; a;, y| A', 3^i, 3^2) 

< Gq,(1, Z2; 1, Z12, Z2\x,y,,y2) 

September 21, 2011 Submitted to IEEE Trans, on Information Theory 



XIA, ZHANG AND DUMITRESCU, ENERGY-EFFICIENT FULL DIVERSITY COLLABORATIVE UNITARY CODE DESIGNS 36 



(1 + a^iEn + E,)){1 + a^{Eu + E^)) ' 

or 



(82d) 



(82e) 



(1 + a^{E, + E2i))(l + + E22)) ■ 

b) If neither Zii,Zi2 nor ^21,^22 belong to the same group, then, by the pigeonhole 
principle, there exists an Xq E X such that one of the following four statements must 
be true: 

i) Zn e ^i.xo and Z21 E 22,xo- As a result, we attain 

Gc,{X,y^,y2) = min G^{x,y;x,y\X,y^,y2) 



< G^iil, Zi, Z21; xo, ^11, Z2)\X, 3^1, 3^2) 

„2 



(1 + + E2i))(l + a^{E2 + Eu)) ' 

ii) Z12 E ^i.xo and Z21 E 22,xo- Similarly, we can arrive at 

Go^{x,y^,y2) < G^{{i,z,,Z2i;xo,z,2,Z2)\x,y^,y2) 

8a^ 



(82f) 



(1 + + E2l))(l + ^2(^2 + Ei2)) ' 

iii) Zn G Zi,^^ and Z22 G 22,xo- Then, 

^^(^,3^1,3^2) < G„((l,Zi,Z22; 0:0,^11,^2)1^, 3^1, 3^2) 

8«2 



(82g) 



(1 + a^{Ei + ^22))(1 + a^{E2 + Eu)) ' 
iv) Z12 E Zi^xo and Z22 E Z2,xo- Also, we can have 

Ga(A',3^i,3^2) < G„((l,Zi,Z22;xo,Zi2,Z2)|A',3^i,3^2) 

8a2 



(82h) 



;i + a^iEi + ^22))(1 + "2(^2 + ^12)) ' 



(82i) 



Now, comparing ([82]) results in a common upper bound on Ga{X, 3^i, 3^2) in Case 2 as follows: 

„2 



G'«(A',3^i,3^2) <max 



to 



(1 + a^iEn + ^2))(1 + "2(^12 + ^2)) ' 

8a2 

(1 + a\Ei + E2i)){l + a^{Ei + E22)) ' 
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6a 



max 





(1 + a^{En + ^2))(1 + a^{Ei2 + ^2)) 
8«2 



;i + a^{E^ + E22))(l + a2(E2 + E12)) ' ' ^^^^ 
since E2 + £"12 < Ei + E21, E12 + E2 < Ei + E22 and E2 + En < Ei + E21 by Lemma |3j 
Following the same trick as the case when 6=1, using the geometrical and arithmetical mean 
inequality first and then, Lemma |4] arrive at the fact that 

8 



;i + a^iEn + E2)){1 + a%Ei2 + E2)) 



(1 + a^{Ei + E22)){1 + a\E2 + £12)) 
Therefore, it follows from this that 



< 



< 



Ef' + eS + + m 



2' 



2 ■ 



GIA", 3^1, 3^2) = max G^A-, 3^1, 3^2) 



< max max < r— ^ ,, , — — — -- , 

- c \{l + a^{En + E2)){l + a^{Ei2 + E2)) 

8a2 ^ 
(1 + a^{E^ + E22)){1 + a'[E2 + £12)) J 

< max 



Ef^ + + JEf^ + Ef^ 



2 ' 



Y, if r is even, 
if r is odd, 



where we have used Lemma [5] in the last step. All the above discussions can be conclude as the 
following property: 
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Tilt T? Q to 




1 y\ 
1^1 
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1 


U.ZJU 






U.J 




U. iZ / 


z 


S OA Ayr 8 OAA/T 


U.ZJH- 


1 S 


V.Vojy 


z 


8 n A A/T 1 A O A Ayr 


U.ZUD 


1 7^^ 




z 


1 A r» A Ayr 1 n A A/T 


xj.Voy 


z 




9 
z 


1 A OA Ayr ^9 OAA/f 


n 1 ^7 

U. 1 J / 


2.25 


0.0294 


2 




0.121 


2.5 


0.0156 


1 


32-QAM, 32-QAM 


0.125 


2.75 


0.0116 


2 


64-QAM, 64-QAM 


0.0762 


3 


0.00820 


2 


64-QAM, 128-QAM 


0.0640 


3.25 


0.00633 


2 


128-QAM, 128-QAM 


0.0563 



TABLE I 

Maximum coding gains for different transmission bit rates using optimal designs 



Property 3: For given r and X = {1, — 1, j, —j}, we have 

^ ^, if r is even, 

(84) 



2 ; 



if r is odd. 



Properties [TJ |2] and [3] lead us to giving tlie following theorem as one of the optimal solutions 
to Problem [21 

Theorem 3: One of the optimal solutions to Problem [2] is given as follows: 

(1) If r = 4, then, 6 = 0,p = q = 2, X = {l},yi = % = Zi = Z2 is the 4-QAM constellation 
and 5 = I/V2. Moreover, the optimal coding gain is Ga{X,yi,y2) = 1/4. 

(2) If r = 10, then, 5 = 0,p = q = 5,X = {l},5^i = 5^2 = = ^2 is the 32-QAM 
constellation and 5 = ^4^ . Moreover, the optimal coding gain is G5(A:', 3^i, 3^2) = 



(V^+V^)2 ■ 

(3) If r is even, then, p = {r + 2)/2,q = r/2, yi = 5^i,opt and 



if r = 6, 



a 



, ^ ^ ^ ^ if r 7^ 6, 10. 



September 21, 2011 



Submitted to IEEE Trans, on Information Theory 



XIA, ZHANG AND DUMITRESCU, ENERGY-EFFICIENT FULL DIVERSITY COLLABORATIVE UNITARY CODE DESIGNS 



39 



Furthermore, the optimal coding gain is given by 



^ 8 if r = 6, 

if r=^ 6,10. 



(86) 



(4) If r is an odd integer exceeding 4, then, 5 = l,p = q= (r + l)/2, = 3^j,opt and 

^ if r = 9, 

if r = 5, 



a 



t/{E[''+Ei'W,''>+Ei\') 



Furthermore, the optimal coding gain is given by 



if r 7^ 5,9. 



8 



if r = 9, 
if r = 5, 



if r ^ 5,9. 



(87) 



(88) 



The proof of Theorem [3] is given in Appendix |Ej The maximum coding gains using the optimal 
UFCPs determined by Theorem [3] are listed in Table |I] for various transmission bit rates, which 
is also shown in Fig. [2} Some observations on Theorem [3] are made as follows: 

(1) Theorem |3] tells us that the training scheme based on the Alamouti code using the 4-QAM 
and 32-QAM constellations is optimal when either one bit or 2.5 bits per channel use is 
transmitted. 

(2) In spite of the fact that from Proposition [5] we know that increasing the number of 
the groups is increasing the minimum Euclidian distance of the constellation y, the 
accumulated minimum Euclidian distance along the two transmitter antennas between two 
distinct groups is always equal to 8. In addition, increasing the number of groups is also 
increasing the size of the constellations and thus, increasing the energies of the three corner 
points. As a result, the UFCP code using four groups cannot enable the optimal coding 
gain. 
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VI. Simulations 



In this section, we carry out computer simulations and compare the error performance of the 
unitary UFCP code design proposed in this paper with those of other schemes in the literatures 
which can be used in a small noncoherent MISO system having two transmitter antennas and 
a single receiver antenna, where channel state information is completely unknown at both the 
transmitter and the receiver and the coherence time is T = 2M = 4. All the schemes that we 
would like to compare here are described as follows: 

(a) Differential unitary code based on Alamouti coding scheme and PSK constellations. This 



design with the fast closed-form ML decoder was proposed in [57 1, [59 1 and two unitary 
codeword matrices are Ui = I2 and 



1 

71 



Sl S2 



X 



(89) 



where si and S2 are randomly, independently and equally likely chosen from the T'^-dxy and 
2'"2-ary phase shift keying (PSK) constellations, respectively, with the two integers a and h 
determined as follows: 

'^^ = ^2 = ^ if r is even, 

2 _ (90) 
ri = ^,^^2 = ^ if r is odd. 

For the necessity of performance comparison and decoding with the GLRT receiver, these two 
unitary matrices are normalized and then, stacked into one codeword matrix, which is denoted 
by S„, 



1 

71 



X 




(91) 



where the normalization constant assures E[tr(S^Sa)] = 2. 

(b) SNR-efficient training Alamouti code. This SNR-efficient training scheme using the 



Alamouti code was presented in [89|. The codeword matrices are characterized by 

/ /E^ \ 







(92) 
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where si and S2 are randomly and equally likely chosen from either the 2'^^-ary and 2^"^-axy 
PSK constellations or cross QAM constellations, respectively, with the determination of the two 



integers ri and r2 being the same as ( [90| ). The energy constant Eb is normalized in such a way 
that E[tr(S^Sfc)] = 2. Here, the optimal average energy distribution over the training phase and 
communication phase is attained by maximizing the training efficiency [[5|, [43 1, |89|. 

(c) Optimal unitary UFCP code. The code design is proposed in this paper and the codeword 
matrix is of the form: 

fx \ 



1 



1 + a^l^/iP + a'^\y2\ 



xeX,yie 3^1,^/2 e 3^2, 



(93) 



X 
ayi ay2 
\ -ayi ayi J 

where the optimal energy scale a and three constellations X, and 3^2 are determined according 
to Theorem [3l 

It can be seen that the above three transmission schemes have the same spectrum efficiency, 
i.e., each transmission rate is = r/4 bits per channel use. To make all error performance 
comparisons fair, we decode all the codes using the GLRT detector, i.e., 

S = arg max Tr (t^S (S^S) S^t) . 

All the average codeword error rates against SNR are shown Fig. [3} It is observed that the 
optimal unitary UFCP code designed in this paper performs the best error performance among 
all the three coding schemes. 



VII. Conclusion And Discussions 

In this paper, we have considered a wireless communication system having two transmitter 
antennas and a single receiver antenna, in which the channel coefficients are assumed to be 
unknown at either the transmitter or the receiver, but remain constant for the first 4 time slots, 
after which they change to new independent values that are fixed for the next 4 time slots, and 
so on. For such a system, we have developed a novel concept called the uniquely factorable 
constellation pair for the systematic design of full diversity unitary space-time block code. By 
simply normalizing the two Alamouti codes and carefully selecting three constellations, a full 
diversity unitary code design with a symbol rate 3/4 has been attained. It has been shown that 
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it is the unique factorization of constellation pairs that guarantees that the unique identification 
of both the channel coefficients and the transmitted signals in the noise-free case as well as 
full diversity in the noise case. In other words, both the unique identification and full diversity 
require that the constellation pair must be designed in such a cooperative way that factorization 
in the product sense is unique-able. It is for this reason that we have named the code proposed 
in this paper as the UFCP code. In addition, to further enhance error performance, the optimal 
unitary UFCP code enabling the maximum coding gain has been designed from a pair of energy- 
efficient cross QAM constellations subject to a bit rate constraint. After a careful examination 
of the fractional coding gain function, in this paper we have taken two major steps maximizing 
the coding gain: 

1) The energy scale has been carefully designed to compress the first three largest energy 
points of the QAM constellations in the denominator of the objective; 

2) The two UFCPs have been designed so carefully that the one constellation collaborates 
with the other two constellations through the two transmitter antennas maximizing the 
minimum of the numerator and at the same time, avoiding the comer points with the 
largest energy as many as possible achieving the minimum. 

In other words, the optimal coding gain has been obtained by constellations collaboration and 
energy compression. It is for this reason that we have also called the optimal UFCP code 
designed in this paper as the energy-efficient collaborative UFCP code. Computer simulations 
have demonstrated that error performance of the optimal unitary UFCP code presented in this 
paper outperforms those of the differential code and the SNR-efficient training code, which, to 
the best knowledge of the authors, is the best code in current literatures for the system. 

As we have seen, the concept of the UFCP plays an important role in the systematic design 
of energy-efficient full diversity unitary space-time block codes for the small MIMO system 
having the two transmitter antennas and a single receiver antenna. However, the constructions 
and properties on the UFCP and the related transmission scheme which have been reported in 
this paper are just initiative. Some significant issues still remain unsolved: 

1) The construction of the optimal UFCPs for the design of the unitary space-time block 
code has been derived from the cross QAM constellations. How about the hexagonal 
constellations? since the hexagonal constellations carved from the Eisenstein integer ring 
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are supposed to be more energy-efficient tlian tlie QAM constellations carved from the 
Gaussian integer ring [j86j . Generally, which constellation is optimal to generate a unitary 
UFCP space-time block code with the optimal coding gain? 

2) Instead of a pair of coprime PSK constellations, whether is the UFCP constructed in this 
paper used to systematically design full diversity noncoherent space-time block codes for 
a general MIMO system by following the way similar to [74]? 

3) The coding scheme which has been adopted in this paper is the Alamouti scheme. In spite 
of the fact that the Alamouti code is optimal in many senses for such coherent system, it is 
not optimal anymore for such noncoherent system, since it was proved that unitary codes 
are optimal for general noncoherent MIMO communications, whereas the Alamouti code 
resulting from the QAM constellation is not unitary in general. Only when the constellation 
is the PSK, the resulting Alamouti code is unitary. However, the PSK constellation is not 
as energy-efficient as the QAM constellation. In addition, although this paper has proposed 
a simple method for the design of the unitary code jusy by normalizing the two Alamouti 
codes, a deep insight into the fractional coding gain function exposures the drawback of 
the Alamouti scheme in the nocoherent case, i.e.. Too large energies are contributed to the 
denominator. Hence, a question is: is it possible to find another coding scheme that has 
the same minimum of the numerator but a smaller maximum of the denominator as the 
Alamouti scheme? 

This paper has just casted a brick so that the jade may be attracted. 

Appendix 

A. Proof of Proposition |5] 

We consider the following two cases: 

Case 1: \X\ = 2, i.e., there are only two elements in X. Let X — {xi,X2}, where Xi,X2 G 
{1, —1, — By Proposition |5} we have Z = Z^^ U Z^^ with Z^.^ fl Z.^.^ = $. In fact, Z^^ = 
xiy = x^y and thus, dmmiZ,j;J = d^iniZ^^) = (imin(3^)- Let P be one of the corner point 
in Z with the largest energy. Without loss of generality, we can always assume P E Z^^. The 
following three possibilities need to be considered seperately. 

1) K = ?). In this case, let P\ is the first nearest neighbor of P, i.e., \P — P\\ =2, and P2 
is the second nearest neighbor of P, i.e, |P — P2I = 2a/2. If Pi G Z.^.^, then, (imin(3^) = 
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rfmin(2..J <\P- Pi\ < 2; If P2 e then, dr^Uy) = rfmin(^.J < \P ' ^2! < 2v^; 
Otherwise, both Pi and P2 must lie in Z^^ and as a result, dmmiy) = d^miZx^) < 
\Pi — -P2I = 2v^. At any rate, the minimum distance of 3^ is upper-bounded by 

d^n{y) < (94) 

On the other hand, it can be verified directly by calculation that the constellation given 
by ^ satisfies three conditions: (a) |3^S| = 4; (b) rfmi„(3^ipl) = 2v^; (c) and X^^, 
indeed form a UFCP. Therefore, J^^pi is optimal. 
2) = 5. In this case, P has the two nearest neighbors, which are denoted by Pi and P2, 
i.e., \P — Pi\ = |P — P2I = 2. If one of Pi and P2, say. Pi, belongs to the Group-xi, Z^^^, 
then, d^iniy) = d^miZ^-^) < |P — Pi| =2. Otherwise, both Pi and P2 must lie in Z^^ 
and as a result, d^nmiy) = '^min(^x2) ^ l-Pi — -P2I = 2\/2. Hence, the minimum distance 
of 3^ is always upper-bounded by 

dmUy) < 2V2. (95) 

Now, following the argument similar to the Possibility 1) of = 3, we can say that the 
constellation y^lj^ given by ([T5]) is indeed optimal. 



opt 

3) K > 4 is even. Similar to the possibility of K = 5, we can prove d^ui{y) < 2\/2. In 
addition, notice that the constellation determined by ([16]) has two properties: (a) {y^ptl = 
2^; (b) yi—y2 = 2(1 — for any two distinct points yi and y2 in 3^opt' where z is some 
complex integer. As a consequence, (imin(3^) = 2v^. Let us now check whether or not 
such a pair of X^^^ and y'^^^_ constitutes a UFCP. Suppose that there exist x,x ^ X and 
y,y E y such that xy = xy, then, (y — y)x = {x — x)y. If x = x, then, y = y, since x 0. 
\f X ^x, then, x — x = ±(1 — j). Combining this with the fact {y — y) = 2(1 — j)z, we 
have 2z = ±y, which is impossible. Thus, x = x and y = y. In other words, a pair of the 
constellations X^ipl and y^^l constitutes a UFCP. Therefore, in this case, yj^^^ is optimal. 

4) K > 5 is odd. Following almost the same discussion as the possibility of the even K 



exceeding 2, we can also arrive at the fact that the constellation given by (Tf) is still 
optimal. 

Case 2: \X\ = 4. In this case, possibilities for K = 3 and 5 can be verified directly by calculation. 
In addition, since the possibility for even K greater than 2 is similar to that for odd K exceeding 
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5, here we only provide a proof for the situation when K is an odd number greater than 5. By 
Proposition Q Z = ZiU Z^iU Zj U Z^j with n Z^., = $ for A; 7^ £ = 1, 2, 3, 4. Actually, 
= ^ly^ since Xk E X = {±l,±j}. Thus, we have dmmiZ^t) = c?mm(3^)- Let Pi for 
i = 1, - ■ ■ , 9 be nine points in the first quadrant around the comer of Z shown in Fig. |4j Using 
the pigeonhole principle, there exists one Goup, say, Z^,,^ , including at least three of these nine 
points, say, Pi^,Pi^ and Pi.^. Among these three points, if there exist two of them lying either in 
the same row or in the same column, then, d^nn{y) = '^min(^xfe(,) < 4. Otherwise, these three 
points locate in different rows and different columns and thus, d^i^(y) = dminiZxi^^J < 2\/2. 
Therefore, in any case, we can always have djamiy) < 4. On the other hand, notice that the 
constellation determined by ([19]) possesses two features: (a) |3^optl = 2^; (b) i/i — 1/2 = for 
any two distinct points yi and ?/2 in 3^opt' where z is some complex integer. Hence, we have 
dmin{y) = 4. Let us now examine whether such a pair of the constellations X^J^ and forms 
a UFCP. Suppose that there exist x, x G A* and y,y E y such that xy = xy, then, we have 

{y - y)x = (x - x)y. (96) 

If X = x, then, y = y, since x ^ 0. If x ^ x, then, x — x = ±2, ±2j, ±1 ± j. No matter whatever 
situation occurs, once we have substituted itj—y) = Az into ( |96l ), we can always obtain y = 2zo 
for some complex integer zq, where we have used the fact that 2 = (1 + j)(l — j)- This implies 
that y is an even number, which is impossible. Thus, x = x and y = y and X^^l and yj;^^ 
indeed constitute a UFCP. Therefore, in this case, y^^^ is optimal. This completes the proof of 
Proposition [5} □ 

B. Proof of Lemma |3] 

Since the proof of ( |35b| ) is much similar to that of ( |35a| ), we only provide a proof for ([35a]), 
which can be fulfilled by considering the following possibilities: 



1) f = 3. In this case, Lemma[2]tells us that E2 = 10 and E21 = 2. Hence, the inequality ([35a)) 
is reduced to 

8 + Eu<Ei. (97) 



If u = 3, El = 10, Ell = 2, so El - En = 8, inequality ([97]) is hold. If u is an odd 
number exceeding 3, then, by Lemma pi we have Ei = (2^ — 1)^ + (3 x 2^ — 1)^ 
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and Eu = (2 V - 3)^ + (3 x 2 V - 1)2. Thus, we obtain - En = (2V - 1)2 _ 
(2^ — 3)^ = 4 X 2^ — 8 > 8, since 2^ > 4. Hence, in this case, the inequality ( [97] ) is 



also true. If u is an even integer, then, by Lemma |2j again, we have Ei = 2(2 2 — 1)2 and 

En = (2t - 1)2 + (2t - 3)2 so that E^ - En = (2^ - 1)2 _ (2t - 3)^ = 4 x 2t - 8 > 8, 



since u > 4. Hence, the inequality (97) is still true. 



2) V is an odd number exceeding 3. In this case. Lemma [2] gives us that E2 = (2 2^ — 

If + (3x2^- 1)2 and E21 = (2^^ - 3)^ + (3 x 2^^ - 1)^ so that E^ - E^i = 



(2V - 1)2 _ (2V - 3)2 = 4 X (2V - 2). On the other hand, if u is an odd number 
not less than v, then. Lemma [2] gives us that Ei = (2^ — 1)^ + (3 x 2^ — 1)^ and 
En = (2^^ - 3)2 + (3 X 2"^ - 1)^ so that E, - En = (2^^ - 1)^ - (2^^ - 3)^ = 
4 X (2^ — 2) > 4 X (2^ — 2), since the exponential function 2* is increasing and 
u > V. Therefore, in this case, we have Ei — En > E2 — E21, which is equivalent to 
the fact that E2 + En < Ei + E21. If u is an even number not less than v, then, using 
Lemma g again yields Ei = 2(2t - 1)2 and En = (2t - 1)2 + (2t - 3)^. Hence, we have 
El - En = (2^ - 1)2 - (2t - 3)2 = 4 X (2t - 2) > 4 X (2^^ - 2) = E2 - E21. Hence, 
in this case, the inequality p5a[ ) holds. 



3) V is an even number greater than or equal to 2. Similarly, by Lemma |2} we can have 

- E21 = 4 X (25 - 2) and 

{4 X (2t — 2) > E2 — E21 if u is even and not less than 2, 
8 if M = 3, 

4 X (2^^ — 2) > i?2 — £"21 if u is an odd number exceeding 3. 
Notice that when m = 3, f must be 2, since u > v > 2. Then, in this case, E2 — E21 = 
2 — 2 = 0. Therefore, we can always obtain Ei — En > E2 — E21. 
This completes the proof of Lemma |3j □ 

C. Proof of Lemma |4] 



Here, we only give the proof of P6b[ ), since the proofs of the other inequalities are much 
similar. To do that, let us consider the following situations. 

1) w is an even number not less than 4. In this case, since u + v = w, both u and v are either 



even or odd. If both u and v are even, then. Lemma [2] provides us with Ei =2(22 —\' 



2 
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and E21 = E22 = (2t - l)^ + (2! - 3^. Then, we have Ei + E21 = 2(2t - 1)2 + (2f - 
1)2 + (2f - 3)2 = 2(2t - 1)2 + (2^ - 1)2 + (2^ - 3)2. Since u>v,w>u> w/2. 
Now, consider the following function in terms of variable t: 



fit) = 2(25 _ 1)2 + (2^ _ 1)2 + (2^ - 3)2 for t > w/2. 
Then, the first order derivative of f(t) with respect to t is given by 

f'{t) = 2(22 - 1)22 ln2 - (2— - 1)2" ln2 - (2" - 3)2" ln2 for w > t > w/2. 



w — t 



Since the exponential function 22 is increasing and t > w — t, (22 — 1)22 > (2 2 —3)2 
and thus, f'(t) > 0, showing that f{t) is an increasing function. Therefore, we have 



fip) > /(^/2), i.e.. El + E21 > El + E21. Consequently, the inequality (36b) holds in 



this case. Similarly, if both p and q are odd, we can also prove that P6b[ ) still holds. 
2) w is an odd number exceeding 4. Since u + v = w, either u is even and v is odd or u is 
odd and v is even. If u is even and v is odd, by Lemma [2} we obtain Ei = 2(2t — 1)2 
E21 = (2^^-3)2 + (3x2^^-1)2 and hence, + = 2(2t - 1)2 + (2^2^ -3)2 + (3 x 
2^^ - 1)2 = 2(25 - 1)2 + (2^ - 3)2 + (3 X 2"^ - 1)2, where w > u > (w + l)/2. 
This leads us to considering a function: 

g{t) = 2(25 - 1)2 + (2^^ - 3)2 + (3 X 2^^ - 1)2 forw>t>{w + l)/2. 
The first order derivative of g{t) is 

, t ^ t ,u! — t — l ,w — t — l^ , vj — t — Z — £ — 3, 

git) = 2(25 _ 1)25 In2 - (2^^ - 3)2^^ ln2 - 3 x (3 x 2^^ - 1)2^^ ln2 

for w > t > (w + l)/2. Since 25 is increasing and t > w — t, we arrive at the fact that 
I X (25 - 1)25 > I X (2 V - 1)2"^ > (2^^ - 3)2^^. In addition, notice that 
(2 - 1) X (2i - l)2l > (2 - I) X (2^ - 1)2^ = 3v^ x (2v^ x 2^ - 1)2^ > 
3 X (3 X 2^^^^ - 1)2^^^^ because of the fact that 2^^^^ > 1. As a result, g'{t) > 
and git) is increasing. Hence, giu) > gHw + l)/2), i.e., i?i + £'21 > Ei + E21. This 
can be also proved to be true if u is odd and v is even. Analogously, we can prove that 



El + E22 > Ei + E22- Therefore, the inequality p6b[ ) holds. 
This completes the proof of Lemma |4j □ 
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D. Proof of Lemma |5] 

We prove Lemma [5] by considering the following four different cases for w. 

I) w = M, where £ is a positive integer not less than 1. In this case, by the definition of u 



and V in p4\ , we have u = v = w/2 = 2i and thus, Ei = E2, En = E22, which implies 
that El — Ell — E2 — E22- 
2) w = U + 1. Then, (|34j) tells us that n = (w + l)/2 = 2£ + 1 is odd, whereas v = 
(if — l)/2 = 2£ is even. Particularly when i = 1, u = 3 and -y = 2. In this specular case, 
El-Ell = 10-2 = 6 and E2-E22 = 2-2 = 0, resulting in Ei-En > E2-E22. If £ > 1, 
then, byLemma|2|we have Ei-En = 4x 2^^ -8 and ^2-^22 = (2i -l)2-(2i -3)^ = 
4x2^-8. Since ^ = | = £, in this case we obtain Ei - En = E2 - -E22. 



3) w = 4£ + 2. From the definition of u and v given in we know that u = v = 
w/2 = 21+1. This means that both u and v are odd. Then, Lemma [2] gives us Ei — En = 
(2^_l)2_(2^-3)2 = 4x2^-8 and ^2-^22 = (3 x 2^^ - 1)^ - (3 x 2^ -3)2 = 
12 X 2^ - 8. Since 12 x 2^ = 6 x 2^ > 4 x 2^, implying Ei - En < E2 - E22. 

4) w = A£ + 3. In this case, u = {w + l)/2 = 2£ + 2 is even and v = {w - l)/2 = 2i + 1 is 
odd. Specially for £ = 1, we have m = 4, -u = 3 and thus, Ei — En = 18 — 10 = 8 and 
E2—E22 = 10— 2 = 8. Hence, we have Ei—En = E2—E22 in this particular case. If £ > 1, 
then, by Lemma[2]again, we can attain that i^i— Ell = (2^— 1)^- (2t— 3)^ =4x2^— 8 and 
E2-E22 = 12x2^-8. Because ofthe fact that 4x2t = 4x2^ = 16x2^ > 12x2^, 
we have Ei — En > E2 — -^22- 

Now, summing up all the above results ends the proof of Lemma [5| □ 

E. Proof of Theorem |i] 

We prove Theorem [3] by considering the following situations: 
1) r = 4. In this case, we know from Properties [T] [2] and |3] that 

1 

4' 
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Therefore, the optimal coding gain is 1/4, 5 = 0,p = g = 2, and X and y are the 4-QAM 
constellation. 

2) r = 10. Similarly, Properties [1} [2] and |3] tell us that 

Gix^'\y['\yi'^) = 



00 

G{x^-\yf\m < 1. 

Hence, the optimal coding gain is ^ ^ S = 0, p = g = 5, and X and y are the cross 
32-QAM constellation. 

3) r is an even integer exceeding 4 and not equal to 10. In this case, we consider the following 
two possibilities: 
a) If r = 6, then. Properties [T] |2] and [3] give us 

Gix^'\yi'\yi'^) =- ^ 



{V3 + V5) 



2 ' 



Gix^'\yi'\yi'^) =- 



2 ' 



2- 



8 



(75 + ^7) 

Thus, in this case, we have S = l,p = 4, g = 3. 
b) If r 7^ 6, then, we have from Properties [T] [2] and |3] that 

Recall the following facts: 

~(2) ~(2) ~(2) ~(2) ~(2) 

i) E{ ' = E2 ,E{2 = , El ' is the largest energy among all the points in the 
2''/2+i.ary cross QAM constellation and e[^ is the energy of the second neighbor 
of the point with the largest energy in the constellation. 



E['^ + eS + JEi'^ + Eg) 



2 ' 



2 ■ 
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ii) e[^^ is the largest energy among all the points in the 2'/^+^-ary cross QAM 
constellation, while e!^^^ and £'22^ the respective energies of the first and 
second neighbors of the point with the largest energy in the 2''/^-ary cross QAM 
constellation. 

iii) Ef'^ = E^\ E^'^ is the largest energy among all the points in the 2^/^-ary cross 
QAM constellation and Efi is the energy of the first neighbor of the point with 
the largest energy in the constellation. 

Hence, we obtain ^jj > E^^^ > E^l^ and as a resuk, G{X'^^\y[^\y^2^^) > 
G{X^^\y[^\yP). In addition, since 



G{x('\yi'\yi'^) ( \Ie? + E^^ + + E,, 



G{xw, y[\ yn \ J2e[') + + V2Ef ) + 2M?) 



< 1, (98) 



21 



where we have used the facts that E^°^ = > e!^^^ and that 2E[°^ > ^ if 

r > 6, we attain G{X^^\y[^\yi^^) > \ J^f )). Therefore, in this case, 

(1) 



the optimal coding gain is G{X,yi,y2) is G{X^^\yl'' ,y. 



2 



4) r is an odd integer exceeding 4. There are three cases which need to be considered: a) 
r = 5, b) r = 9 and c) r 7^ 5, 9. However, since the discussions on Cases a) and b) are 
much similar to the previous cases r = 6 and r = 10, we only consider Case c). In this 
situation. Properties [T| |2] and |3] provide us with 



G{x('\y^f\yS'^) = 
G{x^'\yi'\y^'^) = 
G{x^'\yi'\W) < 



4 





+ + y 

8 


/£<»> 






8 







Ef^ + + JEf^ + 



2 ■ 



Notice the following facts: 
a) E^^ is the largest energy among all the points in the 2'^''+^^/^-ary cross QAM 
constellation, E\y and E\2 are the energies of the first and second neighbors of 
the points with the largest energies in the 2('+'^)/^-ary cross QAM constellations. 
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b) e[^^ = e[^\e[^^ = _E2i^ is the largest energy among all the points in the 
2(^+i)/2.ary cross QAM constellation and e[\^ is the energy of the first neighbor of 
the point with the largest energy in the constellation. 

c) e[^^ is the largest energy among all the points in the 2'^^+^)/^-ary cross QAM 
constellation, E^^ is the largest energy among all the points in the 2*^''~^)/^-ary cross 
QAM constellation and £^21^ is the energy of the first neighbor of the point with the 
largest energy in the 2(''^^)/^-ary cross QAM constellation. 

Now, following the way similar to the previous case where r > 8 is even but not equal 
to 10, we can prove that G{X'-^\y[^\y2 '') is still the optimal coding gain. 
This completes the proof of Theorem [3| □ 
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(a) Diagonal line search over the square QAM (b) Diagonal search over the cross QAM constella- 
constellation for the optimal 3^opi: O tion for the optimal y^^^: O 
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(c) Horizontal and vertical lines search over the (d) Horizontal and vertical lines search over the 
square QAM constellation for the optimal J^^p^: cross QAM constellation for the optimal y^^^: O 

Fig. 1 

The optimal UFCP selections from the square 64QAM and cross 128QAM constellations in Proposition[5] 
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(c) Rb = 1.75 bits per channel use (d) Rb = 2.25 bits per channel use 

Fig. 3 

Performance comparison of unitary UFCP codes with currently available noncoherent codes 
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Fig. 4 

Nine corner points in the cross QAM constellation 
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